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PREFACE 


This  report  documents  work  performed  during  the  period  January  1976 
through  April  1977  by  the  University  of  New  Mexico  under  Contract  F29601- 
76-C-0015  with  DET  1 (CEEDO)  ADTC,  Air  Force  Systems  Command,  Tyndall 
Air  Force  Base,  Florida  32403.  Captain  Robert  F.  Olfenbuttel  managed  the 
program. 

This  volume,  which  documents  program  PHASE4,  is  the  fourth  of  four 
volumes  constituting  the  Air  Force  refuse-collection-scheduling  program 
description.  Except  for  Shell's  sorting  algorithm,  used  in  subroutine 
SHLSRT,  all  algorithms  used  in  program  PHASE4  were  developed  and  coded 
by  Harold  J.  Iuzzolino. 


The  report  has  been  reviewed  by  the  Information  Officer  and  is  re- 
leasable to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  available  to  the  general  public,  including  foreign  nations. 
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1 .  BACKGROUND 

In  designing  the  Air  Force  Refuse-Collection  Scheduling  Program  (RCSP), 
the  fundamental  objective  was  to  reduce  collection  costs.  The  most  effective 
way  to  significantly  reduce  cost  is  to  reduce  the  number  of  trips  used  to  ser- 
vice a given  region.  If  a collection  crew  can  be  dropped  from  the  fleet,  the 
cost  of  manpower  will  be  cut.  In  addition,  fuel  and  maintenance  costs  will  be 
lowered  if  the  total  mileage  traveled  by  the  fleet  can  be  reduced. 


2.  OBJECTIVE 

The  primary  purpose  of  program  PHASE4  is  to  produce  maps  and  printed 
schedules  that  allow  for  easy  implementation  of  the  routes  generated  by  the 
RCSP.  PHASE4  also  selects  trips  from  the  two  choices  per  section  produced  by 
PHASE3  in  a manner  that  provides  the  shortest  total  travel  distance.  PHASE4 
then  selects  two  trips  for  each  vehicle  in  a manner  that  gives  the  shortest 
maximum  route  time. 

3.  SCOPE 

This  report  describes  the  workings  of  program  PHASE4.  A program  overview 
is  given,  followed  by  a thorough  description  of  the  logic  involved  in  schedule 
generation.  Flow  charts  provide  a skeleton  of  the  logic  flow.  Input  and  out- 
put files  are  described.  Program  requirements  and  restrictions,  error  messages 
and  error  handling  techniques,  definitions  of  important  variables,  and  an  esti- 
mate of  running  time  are  also  provided. 
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SECTION  II 
PROGRAM  OVERVIEW 

Program  PHASE4  performs  four  basic  tasks.  It  reads  the  input  data  and  per- 
forms some  validity  checking,  pairs  trips  to  minimize  both  the  total  travel  dis- 
tance and  the  maximum  route  time,  prints  the  collection  schedule,  and  produces 
maps  of  the  collection  trips.  When  trip  information  that  does  not  come  from 
program  PHASE3  (e.g.,  data  describing  routes  already  in  use)  is  used  as  input, 

the  two  choices  of  trips  for  each  section,  which  PHASE3  provides,  will  not  be 

■> 

available.  In  this  case,  PHASE4  uses  the  value  of  a trip-pairing  indicator  to 
determine  whether  the  trips  should  be  scheduled  as  single-trip  routes  or  should 
be  paired  sequentially  into  two- trip  routes. 

The  input  to  PHASE4  consists  of  data  files  passed  on  from  previous  pro- 
grams and  data  read  in  from  cards.  The  file  assignments  are  as  follows: 

TAPE1  is  the  segment  data,  TAPE2  is  the  node  data,  TAPE3  is  the  street-name 
data,  and  TAPE9  is  the  path  data.  The  data  to  be  read  from  cards  consist  of  a 
title,  the  unit  of  refuse,  a trip-pairing  indicator,  street  numbers  and  names, 
time  limits  for  trips  and  the  times  and  durations  of  breaks,  vehicle  descrip- 
tions, and  map  bounds.  The  street  numbers  and  names  are  optional  if  they  are 
on  file  TAPE3.  The  map-bounds  record  is  optional  and  may  be  used  to  describe 
which  regions  should  be  plotted  for  a particular  section  and  trip.  If  the 
route  description  is  not  found  on  TAPE9,  then  it  must  be  read  from  cards. 

The  program  consists  of  a main  program  named  PHASE4,  11  subroutines,  and 
2 function  subprograms.  PHASE4  reads  the  title  and  refuse  units  from  cards  and 
the  street  segment  and  node  data  from  files  TAPE1  and  TAPE2.  It  then  calls 
STRIN,  which  looks  for  street  names  on  file  TAPE3.  If  the  street  names  are  not 
found  on  TAPE3,  they  are  read  from  cards.  Control  returns  to  PHASE4,  which 
reads  the  starting  times  and  durations  for  lunch  and  two  breaks  and  the  vehicle 
capacities  and  identifications.  PHASE4  then  looks  for  map-bounds  cards.  If  the 
map  bounds  are  omitted,  travel  will  be  shown  in  only  the  collection  region. 
PHASE4  reads  the  route  description  from  file  TAPE9,  which  has  been  generated  by 
PHASE3. 
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If  route  descriptions  are  not  found  on  TAPE9  or  on  cards,  the  program 
terminates.  Subroutine  CUMTD  is  called  to  determine  the  cumulative  time  and 
distance  for  travel  from  the  garage  to  the  collection  region,  within  the  col- 
lection region,  and  from  the  region  to  the  landfill.  CUMTD  also  totals  the 
refuse  accumulated  for  each  trip.  A summary  of  the  time,  distance,  and  refuse 
quantity  for  each  section  is  printed  by  PHASE4.  The  trips  are  then  ordered  by 
SHLSRT  according  to  vehicle  capacity.  Appropriate  morning  and  afternoon  trips 
are  paired  in  such  a way  as  to  minimize  the  total  distance  traveled.  The  sec- 
tion numbers  comprising  the  routes  are  stored  by  program  PHASE4.  TAPE9  is  then 
rewound,  and  P0SIT9  is  called  to  position  the  tape  at  the  correct  trip  data. 

The  route  data  are  read  from  TAPE9,  and  PRSCHED  is  called  to  print  the  schedule. 
PHASE4  then  scans  the  map-bounds  data.  Each  time  bounds  other  than  zero  are 
found,  PATHPLT  is  called  to  plot  the  appropriate  vehicle  path.  Subroutine 
STNAME  appends  the  street  names  before  the  map  is  drawn.  MAPPLT  then  plots  the 
map,  and  control  returns  to  PHASE4. 

The  flow  of  control  from  one  subprogram  to  another  is  shown  in  Figure  1. 
Within  each  subprogram,  only  the  first  call  to  each  other  subprogram  is  shown. 
(Subroutines  PLOTS,  PLOT,  and  SYMBOL  are  system  subroutines  and  are  not  in- 
cluded in  the  description  of  PHASE4.) 


SECTION  III 
PROGRAM  LOGIC 


The  logic  for  program  PHASE4  is  described  from  three  viewpoints.  The  first 
description  is  task-oriented.  The  second  describes  the  storage  and  use  of  input 
data.  The  third  describes  each  subroutine  and  the  main  program  in  terms  of  its 
purpose  and  the  manipulations  performed  within  it. 


1 . PROGRAM  TASKS 

The  processing  performed  by  program  PHASE4  can  be  grouped  into  four  tasks: 
checking  the  input  data,  pairing  the  trips,  generating  the  printed  schedule,  and 
plotting  the  route  maps. 

Trip  pairing  is  accomplished  by  the  statements  in  the  main  program,  from 
the  loop  through  statement  530  through  statement  620.  Vehicle  capacities  for 
the  trips  are  sorted  into  increasing  order;  the  sequence  numbers  of  the  trips 
are  carried  along  during  the  sort.  The  distances  saved  by  using  the  second 
trip  of  the  day  rather  than  the  first  trip  are  computed  and  sorted  into  increas- 
ing order  within  each  vehicle  capacity.  The  sequence  numbers  of  the  trips  are 
carried  along  during  the  sort.  The  trip  starting  at  the  garage  will  be  used 
for  the  half  of  the  sections  that  provides  the  smaller  distance  savings.  The 
trip  starting  at  the  landfill  will  be  used  for  the  remaining  half  of  the  sec- 
tions. This  procedure  causes  half  of  the  sections  to  be  serviced  by  trips 
starting  at  the  garage  and  half  to  be  serviced  by  trips  starting  at  the  land- 
fill  and  causes  the  total  travel  distance  for  all  sections  to  be  the  minimum 
obtainable  from  the  two  choices  of  trips  for  each  section. 

The  lengths  of  all  trips  starting  at  the  garage  for  a given  capacity  vehi- 
cle are  sorted  into  increasing  order.  The  lengths  of  the  trips  starting  at  the 
landfill  are  also  sorted  into  increasing  order.  The  shortest  trip  starting  at 
the  garage  and  the  longest  trip  starting  at  the  landfill  are  assigned  to  one  ve- 
hicle. The  next  shortest  trip  starting  at  the  garage  and  the  next  longest  trip 
starting  at  the  landfill  are  assigned  to  the  second  vehicle.  This  procedure 
is  repeated  until  all  of  the  trips  have  been  assigned  to  vehicles.  Selection 
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of  trips  in  this  manner  provides  a reasonable  balance  in  route  lengths  by  pro- 
ducing the  shortest  possible  maximum  distance  traveled  and  the  longest  possible 
minimum-length  route  for  each  vehicle. 

The  pairing  procedure  is  repeated  for  all  trips  of  each  vehicle  capacity. 
If  an  odd  number  of  sections  are  serviced  by  a vehicle  of  a particular  capa- 
city, the  longest  trip  starting  at  the  garage  stands  alone. 

Subroutine  PRSCHED  is  called  once  for  each  vehicle  to  produce  a printed 
schedule  of  the  route.  Each  segment  in  the  trip  produces  a line  of  printed 
output  in  the  schedule,  unless  a travel  stretch  is  found  on  pieces  of  the 
same  street  having  the  same  speed  limits.  Lines  indicating  lunch  and  break 
times  are  inserted  in  the  schedule  immediately  following  the  description  of 
the  segment  traversed  during  the  lunch  or  break  starting  time.  Unloading 
time  at  the  landfill  is  also  indicated  by  a separate  line  in  the  schedule. 

The  maps  are  plotted  in  three  parts.  The  main  program  calls  subroutine 
PATHPLT  to  plot  the  path  traveled  by  the  collection  vehicle.  Travel  is  indi- 
cated by  a solid  line,  and  collection  is  indicated  by  a dashed  line.  If  col- 
lection is  from  one  side  of  the  street  at  a time,  additional  dashes  are  plotted 
perpendicular  to  the  direction  of  travel.  Subroutine  PATHPLT  counts  the  number 
of  times  each  segment  is  traversed  in  each  direction.  The  travel  path  is  plot- 
ted on  the  correct  side  of  the  street,  and  subsequent  traversals  are  drawn  clo- 
ser to  the  center  of  the  street.  Arrowheads  indicating  the  direction  of  travel 
are  appended  to  the  path  at  approximately  3-inch  intervals. 

The  main  program  calls  subroutine  STNAME  to  append  street  names  to  the 
map.  STNAME  obtains  a street  name  and  number  and  determines  the  length  of 
space  needed  for  the  name.  All  of  the  segments  are  examined,  and  segment  num- 
bers and  lengths  for  segments  of  the  given  street  that  are  within  the  bounds 
of  the  map  are  saved.  The  segment  lengths  are  sorted  into  decreasing  order, 
with  the  segment  numbers  carried  along  during  the  sort.  If  the  street  name 
will  fit  on  the  longest  untraveled  segment,  the  name  is  plotted  witnin  the 
bounds  of  that  segment.  If  the  name  will  not  fit  in  the  longest  untraveled 
segment,  the  longest  segment  is  examined.  If  the  name  will  fit  there,  it  is 
plotted  outside  the  segment.  This  procedure  is  repeated  for  each  street  name. 
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The  main  program  calls  subroutine  MAPPLT  to  plot  the  sides  of  the  streets. 
Each  segment,  in  the  map  description  is  examined.  If  the  segment  is  at  least 
partially  within  the  map  bounds,  subroutine  SHAPCOM  is  called  with  a 0.1-inch 
displacement  from  center  indicated  for  the  side  of  the  street.  Subroutine 
SHAPCOM  adjusts  the  scale  of  the  segment  so  that  the  side  of  the  segment  is 
plotted  V2/10  (0.1414)  inch  shorter  than  the  actual  length  of  the  segment. 

In  this  way,  sides  of  street  segments  that  meet  at  angles  of  90  degrees  or 
more  are  prevented  from  overlapping. 

2.  DATA  STORAGE 


Program  PHASE4  obtains  data  from  five  sources:  card  input  and  files  TAPE1 , 
TAPE2,  TAPE3,  and  TAPE9.  Files  TAPE1 , TAPE2,  and  TAPE3  are  the  segrent,  node, 
and  street-name  data  produced  by  program  RCINPT.  File  TAPE9  is  the  route  data 
produced  by  program  PHASE3.  A plot  file,  TAPE8,  is  generated  by  program  PHASE4. 
The  plot  file  contains  maps  of  the  routes. 

The  first  executable  statement  in  the  main  program  reads  the  problem  title, 
the  refuse  unit,  and  the  consecutive-trip  pairing  indicator  from  the  first  two 
data  cards.  The  title  is  stored  in  blank  COMMON.  The  refuse  unit  is  stored  in 
array  UNITS  in  COMMON  block  U.  Storage  for  the  pairing  indicator  is  local  to 
the  main  program.  These  values  are  not  changed  during  the  remainder  of  the 
program. 


The  second  executable  statement  reads  the  segment  data  from  file  TAPE1 . 
All  of  the  data  except  variable  AVMD  are  stored  in  blank  COMMON.  Storage  for 
variable  AVMD  is  in  COMMON  block  MAPDATA.  These  data  are  not  changed  during 
the  remainder  of  the  program. 


The  third  executable  statement  reads  the  node  data  from  file  TAPE2.  Vari- 
ables NHTOT  and  TOTREF  are  kept  in  storage  local  to  the  main  program.  The  values 
are  cleared  just  before  the  loop  through  statement  620,  and  the  variables  are 
reused  in  that  loop.  The  remaining  data,  variable  KNODES  and  arrays  NODNUM, 

NBS,  XNOD,  and  YNOD,  are  stored  in  COMMON  block  NODDATA  and  are  not  changed 
during  the  remainder  of  the  program. 
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The  main  program  calls  subroutine  STRIN  to  read  the  street-name  data. 

The  street  names  and  numbers  are  buffered  in  from  file  TAPE3,  or  are  read  from 
the  remaining  cards  in  the  first  record  if  TAPE3  is  empty.  The  names  and  num- 
bers are  stored  in  arrays  in  COMMON  block  STREETS  and  are  not  changed  during 
the  remainder  of  the  program. 

After  control  returns  to  the  main  program,  two  cards  of  time-restriction 
data  are  read.  All  of  the  data  are  stored  in  COMMON  block  TIMES.  The  stop 
time  per  house,  stop  time  per  unit  of  refuse,  and  unloading  time  at  the  land- 
fill remain  unchanged  during  the  remainder  of  the  program.  The  maximum  trip 
time  is  changed  to  4 hours  if  the  value  read  is  not  greater  than  zero.  The 
maximum  route  time  is  changed  to  8 hours  if  the  value  read  is  not  greater  than 
zero.  The  starting  time  is  changed  to  8:00  o'clock  if  the  value  read  is  not 
greater  than  zero.  The  vehicle  speed  during  collection  is  changed  to  5 mph  if 
the  value  read  is  not  greater  than  zero.  The  duration  and  starting  times  for 
lunch  and  for  each  of  the  two  breaks  are  not  changed  during  the  remainder  of 
the  program. 

The  main  program  reads  vehicle-identification  cards  in  the  loop  through 
statement  160.  The  values  are  kept  in  storage  local  to  the  main  program  and 
remain  unchanged  through  the  rest  of  the  program. 

The  main  program  reads  map-bounds  cards  in  the  loop  through  statement 
190.  Except  for  array  NTRP,  the  values  remain  unchanged  through  the  rest  of 

the  program.  If  the  trip  number  is  not  greater  than  zero,  it  is  reset  to  1 at 

statement  185.  Map  bounds  are  computed  in  the  main  program  between  statements 
450  and  480  for  any  trip  that  lacks  a map-bounds  card.  All  map-bounds  vari- 
ables are  kept  in  storage  local  to  the  main  program. 

The  part  of  the  route  data  describing  the  path  from  the  landfill  to  the 
garage  is  read  at  statement  270  in  the  main  program.  If  no  data  are  found  on 
file  TAPE9,  the  data  are  read  from  cards  and  are  written  to  TAPE9.  The  remain- 
ing route  data  are  read  shortly  before  statement  400  of  the  main  program.  If 

the  data  come  from  cards,  they  are  written  to  TAPE9.  All  of  the  variables  are 

in  storage  local  to  the  main  program.  The  count  of  segments,  N,  and  the  sec- 
tion number,  NSC,  are  reused  later  in  the  program.  Values  for  vehicle  capacity 


and  vehicle  load,  TRC  and  TRL,  are  updated  each  time  path  data  are  read.  The 
path  length,  node  numbers,  segment  numbers,  and  col lection-or-travel  indicators 
remain  unchanged  through  the  rest  of  the  program. 


• The  route  data  are  reread  from  TAPE9  in  the  main  program  in  the  loop 

through  statement  640.  The  trips  are  not  printed  in  order  of  occurrence  unless 
the  consecutive-trip  pairing  indicator  is  greater  than  zero.  Therefore,  subrou- 
tine P0SIT9  is  used  to  position  TAPE9  at  the  trip  to  be  printed.  If  IPAIR  is 
greater  than  zero,  each  pair  of  consecutive  trips  is  assigned  to  one  vehicle  in 
the  order  in  which  the  trips  occur. 

3.  PURPOSE  AND  PERFORMANCE 

Ip  this  section  the  simpler  subroutines  are  described  first  so  their  work- 
ings will  be  understood  when  they  are  mentioned  again  in  the  descriptions  of 
the  more  complicated  subroutines.  The  main  program  is  described  last.  Logic 
flowcharts  are  given  in  Appendix  A.  Complete  program  listings  are  provided  in 
Appendix  B.  In  Appendix  C,  the  more  important  FORTRAN  variables  are  described 
for  each  subroutine. 

a.  Function  I FIND 

Function  IFIND  uses  a binary  search  to  find  a given  number  in  an 
array  and  assigns  the  subscript  of  the  number  as  the  value  of  IFIND.  If  the 
number  is  not  found,  the  function  sets  the  value  of  IFIND  equal  to  the  negative 
of  the  subscript  at  which  the  number,  to  be  in  numerical  order,  should  be  in- 
serted. (The  array  is  assumed  to  be  in  increasing  order.)  The  function  has 
three  arguments.  Argument  NUM  is  the  number  that  is  sought  in  array  IARRAY. 
Argument  IARRAY  is  the  array  to  be  searched.  Argument  LEN  is  the  length  of 
array  IARRAY. 

Function  IFIND  begins  by  checking  to  determine  whether  LEN  > 0.  If 
LEN  < 0,  the  function  assigns  a value  of  -1  to  IFIND,  and  control  returns  to  the 
calling  program.  The  value  -1  indicates  that  the  number  sought  is  not  in  the 
array  and  would  be  stored  as  the  first  entry  in  the  array.  The  binary  search 
uses  variables  II,  IP,  and  IF  as  pointers.  II  is  the  subscript  of  the  front 
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of  the  region  being  searched,  IP  is  the  subscript  of  the  item  being  compared 
to  the  number  being  sought,  and  IF  is  the  subscript  of  the  last  item  in  the 
region  being  searched.  Variable  II  is  set  to  1 at  statement  5.  Variable  IF 
is  set  to  the  length  of  the  array.  The  pointer,  IP,  is  the  subscript  midway 
between  II  and  IF  and  is  computed  at  statement  10. 

Following  statement  10,  NUM  is  compared  to  IARRAY(IP).  If  NUM  < 
IARRAY(IP) , control  transfers  to  statement  20,  indicating  that  the  number  is 
in  the  front  half  of  the  region  being  searched.  At  statement  20  the  final 
pointer  is  moved  to  the  subscript  preceding  the  point  just  searched.  If  NUM  > 
IARRAY(IP),  control  transfers  to  statement  30,  indicating  that  the  number  being 
sought  follows  the  subscript  just  inspected.  At  statement  30  the  initial 
pointer,  II,  is  set  to  the  present  pointer,  IP,  plus  1.  If  the  number  sought 
is  found  at  IARRAY(IP),  control  transfers  to  statement  50,  where  IFIND  is  set 
equal  to  the  current  pointer.  Control  returns  to  the  calling  program.  Where 
NUM  is  unequal  to  IARRAY(IP),  the  initial  or  final  pointer  is  moved  and  control 
resumes  at  statement  40.  At  statement  40  the  final  pointer  is  compared  to  the 
initial  pointer.  If  IF  > II,  control  transfers  to  statement  10,  where  the 
search  resumes  on  the  appropriate  half  of  the  region  examined  previously.  If 
the  final  pointer  becomes  less  than  the  initial  pointer,  the  number  sought  is 
not  in  the  table.  In  this  case,  control  resumes  following  statement  40,  and 
the  value  of  IFIND  is  set  to  the  negative  of  the  current  pointer.  If  the  num- 
ber at  the  current  pointer  is  less  than  the  number  being  sought,  IFIND  is  set 
to  -(IP  + 1)  so  the  number  can  be  inserted  in  the  appropriate  place.  Control 
then  returns  to  the  calling  program. 

b.  Subroutine  STRIN 

Subroutine  STRIN  reads  the  street  numbers  and  name;.  STRIN  first 
checks  file  TAPE3  for  street  data.  These  data  would  have  been  written  by  pro- 
gram RCINPT  in  the  first  phase  of  data  entry.  If  the  street  data  are  found  on 
file  TAPE3,  they  are  buffered  in  and  stored  in  core.  The  street  numbers, 

NUMSTR,  and  the  street  names,  NAMSTR,  are  counted  and  totaled  in  variable 
NSTRS.  The  total  numbe”  of  street  numbers  and  names  read  from  TAPE3  is  printed. 
If  no  street  data  are  found  on  TAPE3,  the  street  numbers  and  names  are  read  from 
cards  and  stored.  The  number  of  street  numbers  and  names  read  from  cards  is  then 
printed.  Control  returns  to  the  calling  program. 


c.  Function  HM 

Function  HM  changes  times  of  day  that  have  been  given  in  hours  and 
fractions  of  hours  to  a character  form  of  hours  and  minutes  that  will  be  used 
in  the  printout.  Variable  IH  is  the  number  of  whole  hours,  and  variable  IM  is 
the  number  of  minutes.  When  the  number  of  minutes  in  a given  time  of  day  is 
large  enough  to  be  rounded  off  to  60,  IM  is  set  to  0,  and  one  whole  hour  is 
added  to  variable  IH.  An  ENCODE  statement  is  used  to  produce  the  character 
form  of  hours  and  minutes. 

d.  Subroutine  NUMBER 

Subroutine  NUMBER  appends  numbers  to  plotted  output.  Its  purpose  is 
almost  identical  to  that  of  the  standard  Calcomp  number  routine,  the  primary 
difference  being  that  the  last  argument  in  subroutine  NUMBER  gives  an  alpha- 
numeric format  rather  than  an  integer  format  code. 

Subroutine  NUMBER  has  six  arguments.  The  first  two  give  the  coordi- 
nates, in  plotter  inches,  of  the  lower  left  corner  of  the  field.  The  third 
gives  the  height  of  the  digits,  in  inches.  The  fourth  is  the  number  to  be 
plotted.  The  fifth  is  the  angle  at  which  the  number  is  to  be  plotted,  measured 
in  degrees  counterclockwise  from  the  horizontal.  The  last  argument  is  an  alpha- 
numeric format  up  to  10  characters  long,  which  describes  the  appearance  of  the 
plotted  number. 

Array  TEXT  is  used  to  hold  the  character  representation  of  the  number. 
Up  to  30  characters  are  allowed.  The  first  executable  FORTRAN  statement  sets 
this  array  to  three  words  of  blanks.  The  second  statement  moves  the  format  into 
the  second  word  of  array  FORM.  The  first  and  third  words  of  this  array  have 
been  preset  to  a left  and  a right  parenthesis  by  a DATA  statement.  The  ENCODE 
statement  converts  the  number  from  binary  form  in  variable  NUM  to  character 
form  in  array  TEXT,  according  to  format  FORM. 

A character  count,  variable  NC,  is  set  to  30.  The  loop  through  state- 
ment 10  searches  for  the  last  nonblank  character  in  array  TEXT.  Each  time  a 
blank  is  found,  starting  at  the  end  of  the  TEXT  array,  the  character  count  is 
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decremented  by  1.  When  a nonblank  character  is  encountered,  control  trans- 
fers to  statement  20.  Statement  20  calls  the  standard  SYMBOL  subroutine  to 
plot  the  character  representation  of  the  number.  Control  then  returns  to  the 
calling  program. 

e.  Subroutine  CUMTD 

Subroutine  CUMTD  determines  the  total  distance,  time,  number  of 
houses,  and  refuse  quantity  for  a given  trip.  Argument  ISEG  is  the  array  of 
segments  to  be  covered  in  a particular  trip,  CORT  indicates  whether  collection 
or  only  travel  is  required  on  a particular  segment,  NSG  is  the  number  of  seg- 
ments, D1S  is  the  total  distance,  TIM  is  the  total  time,  NHT  is  the  total 
number  of  houses,  RQ  is  the  total  refuse  quantity,  SCOLL  gives  the  vehicle 
speed  in  the  collection  area,  TSTOPH  is  the  stop  time  per  household,  and 
TSTOPR  is  the  stop  time  per  unit  of  refuse. 

After  the  total  number  of  houses,  the  total  distance,  the  total 
time,  and  the  total  refuse  quantity  are  set  to  0,  CUMTD  checks  each  segment 
and  adds  its  length  to  the  distance  accumulated  to  that  point.  CUMTD  next 
checks  the  collection-or-travel  indicator  for  that  segment.  If  collection  is 
required,  three  totals  are  changed:  (1)  TIM  is  increased  by  the  time  spent 
by  the  vehicle  on  that  segment;  the  time  is  computed  as  the  sum  of  the  length 
of  the  segment  divided  by  the  speed  of  the  vehicle  when  collecting  refuse 
plus  the  estimated  total  collection  time  required  for  that  segment.  (2)  NHT 
is  increased  by  the  number  of  houses  on  that  segment.  (3)  RQ  is  increased  by 
the  amount  of  refuse  for  all  of  the  houses  on  that  segment.  If  no  collection 
is  required  on  the  segment,  only  the  total  time  is  changed.  In  this  case, 
the  total  time  is  increased  by  the  length  of  the  segment  divided  by  the 
allowed  speed  of  the  vehicle  for  that  segment.  Control  then  returns  to  the 
calling  program. 

f.  Subroutine  SHLSRT 

Subroutine  SHLSRT  sorts  an  array  into  either  increasing  or  decreas- 
ing order.  SHLSRT  has  four  arguments:  X,  the  array  to  be  sorted;  A,  which 
is  reordered  as  X is  sorted  so  that  each  A always  corresponds  to  the  same  X; 


. 
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NW,  the  number  of  words  to  be  sorted;  and  SGN,  which  indicates  whether  X is 
to  be  sorted  in  increasing  or  decreasing  order. 


The  pointer  separation,  N,  is  set  at  half  the  number  of  words  to  be 
ordered.  The  pointer  end  value,  K,  is  computed  next.  The  initial  and  final 
pointers,  J and  L,  are  set.  The  values  in  arrays  X and  A at  the  inspection 
point  are  saved  in  variables  XT  and  AT.  SHLSRT  then  determines  whether  the 
final  and  the  initial  location  values  are  in  correct  order.  Variable  SGN 
will  have  a value  of  1 if  X is  to  be  sorted  in  increasing  order,  and  a value 
of  -1  if  X is  to  be  sorted  in  decreasing  order.  If  the  values  are  not  in 
proper  order,  the  initial  value  is  moved  to  the  final  location.  The  final 
pointer  is  set  equal  to  the  initial  pointer,  and  the  initial  pointer  is  moved 
up  by  the  spacing  value,  N. 

SHLSRT  then  determines  whether  the  initial  pointer  has  run  off  the 
front  of  the  array.  If  it  has,  the  saved  values,  XT  and  AT,  are  stored  in 
the  location  designated  by  the  final  pointer.  SHLSRT  next  checks  the  pointer 
spacing.  If  the  spacing  has  not  been  reduced  to  1,  it  is  set  equal  to  half 
the  current  value  and  the  entire  process  is  repeated.  If  the  pointer  spacing 
is  equal  to  1,  control  is  returned  to  the  calling  program. 

g.  Subroutine  P0SIT9 

Subroutine  P0SIT9  positions  file  TAPE9,  the  tape  containing  the 
route-description  data,  at  the  beginning  of  the  trip  to  be  processed.  P0SIT9 
has  two  arguments:  LSEQN  is  the  sequence  number  of  the  last  trip  read  from 
file  TAPE9,  and  NSEQN  is  the  sequence  number  of  the  trip  to  be  processed  next. 


Subroutine  P0SIT9  first  determines  the  position,  on  file  TAPE9,  of 
the  next  trip  to  be  processed.  If  the  next  trip  is  the  next  record  on  the 
file,  control  is  returned  to  the  calling  program.  If  the  next  trip  is  ahead 
of  the  last  trip  data  read  from  the  tape,  TAPE9  is  rewound  and  repositioned 
by  means  of  a dummy  read.  If  records  must  be  skipped  in  order  to  reach  the 
next  trip  data  to  be  processed,  a dummy  read  is  used  to  find  the  correct 
position.  Control  is  then  returned  to  the  calling  program. 


h.  Subroutine  SHAPCOM 


Subroutine  SHAPCOM  sets  up  parameters  in  COMMON  block  COPARM  that  de- 
scribe the  geometric  properties  of  a segment.  These  parameters  are  used  by  sub- 
routine COORD  to  produce  the  coordinates  of  points  on  a segment. 

Subroutine  SHAPCOM  has  six  arguments.  Argument  TOTLEN  gives  the  total 
length  of  the  segment,  in  miles.  Argument  AVMD  gives  the  number  of  miles  per 
map  coordinate  unit  (MCU)  on  the  overall  map.  W is  the  displacement  from  the 
center  of  the  segment  to  the  line  drawn.  DI  and  DF  are  displacements  of  the 
starting  and  ending  points  along  the  segment,  in  miles.  Argument  DIR  gives  the 
direction  of  travel.  The  values  of  the  arguments  are  sent  to  subroutine 
SHAPCOM,  and  all  output  values  from  SHAPCOM  are  placed  in  COMMON  block  COPARM. 

In  COMMON  block  COPARM,  variable  SF  indicates  the  shape  of  the  seg- 
ment. XNI  and  XNF  are  the  x-coordinates  of  the  initial  and  final  nodes  of  the 
segment.  YNI  and  YNF  are  the  y-coordinates  of  these  nodes.  SX  and  SY  are  the 
slopes,  in  MCU  per  mile,  in  the  x and  y directions.  RPR  is  the  reciprocal  of 
the  radius  of  curvature  for  circular  segments  and  the  circular  portions  of  S- 
curves.  Cll  and  C12  are  the  position  differences,  in  MCU,  of  the  starting 
point  and  center  of  a circular  arc  or  the  first  half  of  an  S-curve.  XCTR  and 
YCTR  are  the  center  coordinates,  in  MCU,  for  a circular  arc  or  half  an  S-curve. 
BR1  is  the  distance  in  miles  from  the  start  of  a segment  to  some  particular 
point  on  that  segment.  It  is  not  used  for  straight  segments.  For  circular 
segments,  BR1  is  the  total  perimeter.  For  an  S-curve,  BR1  is  the  perimeter  to 
the  midpoint  of  the  S-curve.  For  a rectangular  segment,  BR1  is  the  distance 
to  the  first  bend  in  the  rectangle.  For  an  angle,  BR1  is  the  distance  to  the 
vertex.  BR2  is  defined  only  for  rectangular  segments  and  angles.  For  a rec- 
tangular segment,  BR2  is  the  distance  in  miles  from  the  beginning  of  the  seg- 
ment to  the  second  bend.  For  an  angle,  BR2  is  the  length  of  the  second  side. 

A value  of  1 is  assigned  to  variable  NPC.  Arrays  A,  B,  RATIO,  XINT,  and  YINT 
are  used  in  various  combinations  according  to  the  segment  shape  being  processed. 

Subroutine  SHAPCOM  begins  execution  by  assuming  that  the  shape  code 
indicates  a straight  line.  BR1  and  BR2  are  set  to  0.  DX  and  DY,  the  x-  and 
y-components  of  the  vector  from  the  initial  node  to  the  final  node  on  the  seg- 
ment, are  computed.  The  x-  and  y-components  of  the  slope  of  the  vector,  measured 


in  MCU  per  mile,  are  computed  and  stored  in  S X and  SY.  Values  are  assigned  to 
A(l),  B(l),  XINT(l),  YINT(l),  and  RATIO(l).  The  shape  code  is  tested;  if  the 
segment  proves  to  be  a straight  line  or  is  not  to  be  plotted,  the  subroutine 
returns  control  to  the  calling  program.  For  any  other  shape  code,  execution 
continues.  The  angle  of  the  vector  from  the  initial  node  to  the  final  node  is 
computed  as  variable  THETA.  The  length  of  the  vector,  0,  is  computed  in  miles. 

The  shape  code  is  checked  again.  If  it  indicates  a shape  other  than  a circular 

arc  or  an  S-curve,  control  transfers  to  statement  40.  If  it  indicates  a circu- 
lar arc  or  an  S-curve,  variables  XE  and  YE  are  set  to  the  x and  y values  of  the 

final  node  of  the  street  segment.  Other  variables  are  set  as  SHAPCOM  assumes 
a circular  arc  shape  code.  The  shape  code  is  then  tested  again.  If  an  S-curve 
is  indicated,  variables  XE  and  YE  are  reset  to  the  coordinates  of  the  midpoint 
of  the  S-curve.  Break  indicator  BR1  is  reset  to  the  perimeter  from  the  starting 
point  to  the  center  of  the  S-curve.  Variable  DO  is  set  to  half  the  distance 
from  the  starting  point  to  the  stopping  point.  Other  variables  necessary  to 
process  an  S-curve  street  segment  are  set.  If  the  shape  code  indicates  either 
a circular  arc  or  a portion  of  an  S-curve  with  a left  direction,  the  value  of 
variable  SGN  is  replaced  by  -SGN.  Other  necessary  variables  are  set,  and  RPR, 
the  reciprocal  of  the  radius  of  curvature,  is  calculated.  The  approximate  RPR 
is  improved  by  a series  of  linear  interpolations.  When  RPR  is  within  the  de- 
sired range  of  accuracy,  the  radius  of  curvature,  R,  is  computed.  A temporary 
variable,  ARG,  is  evaluated.  The  height  of  the  center  of  the  circle  from  the 
line  connecting  the  starting  and  stopping  points  is  set  to  0.  If  variable  ARG 
is  greater  than  0,  the  height,  H,  is  recomputed.  The  distance  to  the  first 
break,  BR1 , is  tested  to  determine  whether  the  circular  arc  is  greater  than  a 
half  circle.  If  it  is,  the  sign  of  the  height  is  changed.  The  x-  and  y- 
coordinates  of  the  center  of  the  circle  are  computed.  The  value  of  RATIO(l) 
is  based  on  the  values  of  W and  RPR.  RATI0(2)  is  set  at  2.0-RATI0(l ) . The 
components  of  the  vector  from  the  center  to  the  starting  point,  Cll  and  Cl 2, 
are  computed.  All  variables  needed  to  compute  points  on  the  circular  arc  or 
the  S-curve  are  now  available,  so  control  returns  to  the  calling  program. 

At  statement  40  the  shape  code  is  tested  to  determine  whether  it  in- 
dicates a rectangular  segment.  For  the  rectangular  segment,  the  distance  from 
the  start  to  the  first  bend,  BR1,  is  computed.  If  this  distance  is  greater 
than  0.05  of  the  total  length  of  the  segment,  control  transfers  to  statement  50, 


where  the  segment  is  treated  as  a rectangular  shape.  Otherwise,  the  rectangle 
is  assumed  to  be  so  shallow  that  a straight-line  approximation  is  adequate,  and 
the  shape  code  is  set  to  0.  Control  then  returns  to  the  calling  program.  At 
statement  50  the  perimeter  to  the  second  bend  in  the  rectangle,  BR2 , is  computed. 
The  values  of  all  other  variables,  including  all  of  arrays  A,  B,  XINT,  YINT, 
and  RATIO,  are  set.  Control  is  then  returned  to  the  calling  program. 

At  this  point,  the  only  segments  that  remain  to  be  processed  are  the 
angles.  Both  BR1 , the  length  of  the  first  side  of  the  angle,  and  BR2,  the 
length  of  the  second  side,  are  computed.  A temporary  variable,  F,  is  computed. 
This  variable  and  the  total  perimeter  of  the  angle  are  used  to  compute  the 
height  of  the  vertex  above  the  line  connecting  the  starting  and  stopping  nodes. 
The  height,  H,  is  then  used  in  the  computation  of  the  x-  and  y-coordi nates  of 
the  vertex  of  the  angle.  All  other  variables  needed  to  compute  points  on  the 
angular  segment  are  computed.  Control  is  returned  to  the  calling  program. 

i.  Subroutine  COORD 

Subroutine  COORD  is  given  a distance,  in  miles,  from  the  beginning 
of  a segment  and  returns  the  coordinates  in  MCU.  Parameters  describing  the 
current  segment  are  stored  in  COMMON  block  COPARM  by  subroutine  SHAPCOM  be- 
fore COORD  is  called.  Argument  CDMLEN  is  the  cumulative  length,  in  miles. 
Arguments  XX  and  YY  are  the  coordinates  returned  for  a point  CDMLEN  miles  from 
the  start  of  the  segment.  IERR  will  have  a value  of  1 or  0,  indicating  whether 
the  values  of  the  coordinates  have  been  calculated  in  subroutine  COORD.  A 
value  of  1 shows  that  an  error  has  occurred  and  that  the  coordinates  have  not 
been  calculated. 

Variable  S is  set  equal  to  the  cumulative  length.  The  shape  code, 

SF,  is  checked.  If  it  is  not  equal  to  0,  control  transfers  to  statement  10. 
Otherwise,  processing  continues  for  a straight  segment.  The  coordinates  of  the 
point  on  a straight-line  segment  are  computed  and  stored  in  variables  XX  and  YY. 
Control  transfers  to  statement  80. 

At  statement  10,  if  the  segment  is  not  straight,  the  shape  code  is 
tested  to  determine  whether  the  segment  is  a circular  arc  or  an  S-curve.  If 
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it  is  neither,  control  transfers  to  statement  30.  For  circular  and  S-curve 
segments,  the  reciprocal  of  the  radius  of  curvature  is  stored  in  RIP.  The 
coordinates  of  the  center  of  the  circular  portion  are  stored  in  XC  and  YC. 

The  components  of  the  vector  from  the  center  of  the  circle  to  the  initial  node 
are  stored  in  Cl  and  C2.  If  the  point  on  the  segment  is  less  than  or  equal  to 
0.999  of  the  first  break  distance  or  if  the  shape  code  indicates  a circular  seg- 
ment, control  transfers  to  statement  20.  The  statements  following  this  test 
change  parameters  to  generate  coordinates  for  the  second  circular  portion  of  an 
S-curve.  The  sign  of  the  reciprocal  of  the  radius  of  curvature  is  then  reversed. 
The  cumulative  distance,  S,  is  set  to  the  distance  from  the  midpoint  of  the  S- 
curve.  The  coordinates  of  the  center  of  the  second  portion  of  the  S-curve,  XC 
and  YC,  are  computed.  Variables  Cl  and  C2  are  recomputed  for  the  new  center. 

The  sine  and  cosine  of  the  angle  subtended  by  the  perimeter  corresponding  to 
S are  computed.  The  coordinates  of  the  point,  XX  and  YY,  are  computed,  and 
control  transfers  to  statement  80. 

At  statement  30,  the  shape  code  is  checked  to  determine  whether  the 
segment  shape  is  a rectangle.  If  it  is  not,  control  transfers  to  statement  60. 
If  it  is,  variable  SGN  is  set  to  1 . If  the  shape  code  indicates  a left  rec- 
tangle, SGN  is  reset  to  -1.  If  S,  the  distance  along  the  rectangle,  is  greater 
than  1.05  times  the  length  of  the  first  side,  control  transfers  to  statement  40 
in  subroutine  COORD.  If  S is  greater  than  0.95  times  the  length  of  the  first 
leg  of  the  rectangle,  S is  set  equal  to  the  length  of  the  first  leg.  The  x- 
and  y-coordi nates  of  the  point  on  the  first  leg  are  computed  by  linear  inter- 
polation, and  control  transfers  to  statement  80. 

At  statement  40,  S is  tested  to  determine  whether  it  falls  on  the 
second  leg  of  the  rectangle.  If  S is  greater  than  1.05  times  BR2,  control 
transfers  to  statement  50.  If  S is  greater  than  0.95  times  8R2,  S is  set 
equal  to  BR2.  The  x-  and  y-coordinates  of  the  point  on  the  second  leg  are 
computed  by  linear  interpolation,  and  control  transfers  to  statement  80. 

At  statement  50,  the  coordinates  of  a point  on  the  third  side  of  the 
rectangle  are  computed.  Control  transfers  to  statement  80. 
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At  statement  60  the  distance,  S,  is  compared  to  the  length  of  the 
first  side  of  an  angle  segment.  If  S is  greater  than  this  length,  control 
transfers  to  statement  70.  If  not,  the  x-  and  y-coordi nates  are  computed  by 
linear  interpolation  for  a point  on  the  first  leg.  Control  transfers  to  state- 
ment 80.  At  statement  70  the  distance  along  the  angle  is  decreased  by  the 
length  of  the  first  leg  of  the  angle.  The  coordinates  of  the  point  on  the 
second  leg  are  computed  by  linear  interpolation. 

At  statement  80,  a linear  interpolation  is  performed  to  obtain,  from 
the  point  on  the  segment,  the  coordinates  of  a point  some  distance  from  the 
center  of  the  street.  The  distance  is  argument  W to  subroutine  SHAPC0M.  Con- 
trol returns  to  the  calling  program. 

j.  Subroutine  PRSCHED 

Subroutine  PRSCHED  examines  the  trip  data  and  prints  the  schedule  for 
a particular  trip.  Argument  NSC  is  the  section  number;  ITRIP,  the  trip  number; 
NT PS,  the  total  number  of  trips  per  vehicle  per  day;  NSP  and  NNP,  arrays  con- 
taining the  numbers  of  the  segments  and  nodes  in  the  path;  C0RT,  the  collection- 
or-travel  indicator;  NSG,  the  number  of  segments  on  a particular  trip;  and  TC, 
the  vehicle  capacity. 

On  the  first  trip  for  each  vehicle,  the  route  printing  is  initialized. 
If  the  current  trip  is  the  first  trip  of  the  day,  PRSCHED  will  print  out  the 
necessary  headings  and  the  time,  in  minutes  and  hours,  at  which  the  truck  leaves 
the  garage.  If  the  current  trip  is  not  the  first  trip  of  the  day,  PRSCHED  will 
print  out  the  time  at  which  the  truck  leaves  the  landfill. 

PRSCHED  then  begins  processing  the  current  segment.  Function  I FIND 
is  used  to  find  the  street  name  and  number.  CORT  is  checked  to  determine 
whether  collection  or  only  travel  is  required  on  that  segment.  According  to 
the  value  of  CORT,  PRSCHED  begins  either  travel  or  collection  processing.  If 
collection  processing  is  indicated,  PRSCHED  determines  whether  there  is  to  be 
pick-up  on  only  the  right  side  or  on  both  sides  of  the  street.  The  totals  for 
the  number  of  houses,  time,  distance,  and  refuse  amount  are  then  accumulated. 

If  only  travel  is  indicated,  the  printed  direction  will  be  to  drive  on  a par- 
ticular street.  Travel  time  and  distance  totals  are  accumulated. 


Subroutine  PRSCHED  next  checks  for  cross  streets  by  using  function 
I FIND  on  the  node  data.  If  there  are  cross  streets,  the  street  is  broken  up 
into  street  segments  by  the  nodes  that  indicate  street  intersections.  The 
direction  to  drive  on  a particular  street,  to  pick  up  on  the  right  side  only, 
or  to  pick  up  on  both  sides  of  the  street,  is  then  printed  out.  At  each  street 
intersection,  the  time  is  checked  to  determine  whether  a morning,  afternoon,  or 
lunch  break  should  have  started  while  the  vehicle  was  either  driving  on  or  col- 
lecting the  last  street  segment.  If  so,  the  break  is  scheduled  for  the  time  at 
which  the  street  segment  is  to  be  completed.  The  direction  to  break  and  the 
beginning  and  ending  times  for  that  particular  break  are  printed  out.  PRSCHED 
next  determines  whether  the  final  street  segment  for  that  trip  has  been  com- 
pleted. If  so,  the  direction  to  unload  is  printed  out,  along  with  the  starting  j 

and  ending  times  for  the  unloading  process.  After  the  last  piece  of  the  trip  j 

has  been  processed,  control  returns  to  the  calling  program. 

i 

k.  Subroutine  STNAME  1 

* 

Subroutine  STNAME  appends  street  names  to  the  maps  of  the  collection 
areas.  STNAME  has  four  arguments:  NSP,  NNP,  CORT,  and  NSG.  NSP  and  NNP  are 
arrays  containing  the  street  segment  numbers  and  the  node  numbers  of  the  travel 
path.  CORT  is  the  collection-or-travel  indicator.  NSG  is  the  number  of  seg- 
ments in  a particular  trip. 

« 

First,  the  map  parameters  are  set.  Then,  a DECODE  statement  in  STNAME 
changes  internally  the  format  of  the  stored  street-name  data.  These  data  are 
examined,  character  by  character.  The  street  names  are  processed  individually 
to  find  the  first  and  last  nonblank  characters.  The  number  of  characters  in 
the  street  name  is  determined  and  stored  in  variable  NCH.  Variable  WIDTH  is 
set  originally  to  half  the  width  of  the  street  on  the  plotted  map.  The  height 
of  the  letters  to  be  printed  is  1.6  times  WIDTH;  therefore,  the  street  name 
will  fit  within  the  street  boundaries.  The  width  of  the  name,  AWDTH,  is  a 
fractional  part  of  the  height  of  the  letters  times  the  number  of  letters. 

A maximum  of  20  segments  with  the  same  street  name  are  saved  tempor- 
arily. The  nodes  defining  these  segments  are  examined  to  determine  whether  the 
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segments  are  within  map  bounds.  If  the  nodes  show  a segment  to  be  within  the 
boundary  of  the  map,  the  segment  number  is  saved  in  array  NSGTEM,  and  the  length 
of  the  segment  is  stored  in  array  FLTEM.  If  the  number  of  segments  saved,  NSV, 
is  less  than  onek  another  street  name  is  checked.  Subroutine  SHLSRT  is  called 
to  sort  the  segments  by  length,  with  the  longest  segment  first.  STNAME  then 
looks  for  an  untraveled  street  segment  long  enough  to  hold  the  street  name.  If 
it  finds  one,  the  name  will  be  printed  inside  the  boundaries  of  that  segment. 
Otherwise,  the  street  name  will  be  written  outside  the  longest  street  segment 
(if  it  is  longer  than  the  street  name). 

The  coordinates  for  printing  the  street  name  are  set  up.  Subroutine 
SHAPCOM  is  called  to  ensure  that  the  printing  of  the  street  name  will  follow 
the  shape  of  the  street  segment.  Next,  the  individual  characters  are  plotted 
by  calls  to  subroutines  COORD  and  SYMBOL.  Subroutine  STNAME  then  determines 
whether  any  more  street  names  must  be  processed.  When  the  processing  of  all 
street  names  has  been  completed,  control  is  returned  to  the  calling  program. 

1.  Subroutine  MAPPLT 


Subroutine  MAPPLT  draws  a street  map  with  double  lines  representing 
the  sides  of  the  streets.  The  subroutine  has  three  arguments.  The  first,  NRT, 
is  the  route  number.  The  second,  ITRIP,  is  the  trip  number.  The  third,  NTPS, 
is  the  number  of  trips  allowed  per  vehicle  per  day. 


The  coordinates  of  the  region  bounding  the  map  are  stored  in  COMMON 
block  MAPDATA.  In  this  COMMON  block,  variables  XMIN  and  XMAX  are  the  minimum 
and  maximum  x-coordinates  for  the  map.  XLEN  is  the  length,  in  inches,  of  the 
map  in  the  x-direction.  YMIN,  YMAX,  and  YLEN  are  the  corresponding  variables 
in  the  y-direction.  YHCUT  is  the  height,  in  plotter  inches,  at  which  the  map 
must  be  sliced  into  strips.  Variable  AVMD  contains  the  miles  per  MCU  conver- 
sion factor  for  each  man.  WIDTH  is  half  the  width  of  the  street,  in  plotter 
inches. 


MAPPLT  begins  by  retrieving  or  computing  the  map  bounds,  the  height 
of  a strip  of  the  map  (PHGT),  the  maximum  length,  the  number  of  map  strips  (MX), 
the  map-scale  factors,  and  the  interval,  in  MCU,  at  which  the  strips  are  to  be 
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The  loop  through  statement  210  controls  the  plotting  of  each  side  of 
the  street.  Variable  W is  the  width,  in  miles,  from  the  center  to  the  side  of 
the  street.  The  loop  through  statement  200  tests  each  segment  to  see  whether 
it  falls  within  the  frame  of  the  map;  if  it  does,  the  segment  will  be  plotted. 
Variables  NI  and  NF  are  set  equal  to  the  numbers  of  the  nodes  bounding  the  seg- 
ment. The  midpoint  coordinates  of  the  segment  are  saved  in  variables  XMD  and 
YMD.  The  lines  in  the  node-number  array  at  which  the  initial  and  final  nodes 
occur  are  saved  in  variables  NS1  and  NS2.  The  coordinates  of  each  node  are  re- 
trieved. 


Initially  the  segment  is  assumed  to  be  entirely  within  the  bounds,  and 
indicators  INBI,  INBM,  and  INBF  are  set  to  1 . If  the  coordinates  of  the  initial 
node  lie  outside  the  frame  of  the  map,  INBI  is  set  to  0.  Similar  tests  are  made 
on  the  coordinates  of  the  midpoint  of  the  segment  ana  the  coordinates  of  the 
final  node  of  the  segment.  If  all  three  points  are  outside  the  frame  of  the 
map,  control  transfers  to  statement  200  and  the  segment  is  not  plotted.  For 
segments  that  are  at  least  partially  within  the  frame  of  the  map,  the  total 
length  of  the  segment,  in  miles,  is  saved  in  variable  T0TLEN.  The  number  of 
points  to  be  used  in  plotting  half  the  segment,  NPMID,  is  computed.  The  number 
will  be  restricted  to  a maximum  of  10  points.  The  total  number  of  points  per 
segment,  NPPSEG,  is  set  to  twice  NPMID. 

Subroutine  SHAPC0M  is  called  to  set  up  the  parameters  needed  to  gen- 
erate coordinates  of  points  on  the  segment.  The  cumulative  length  along  the 
segment  is  initially  set  to  0.  A step  size,  DS,  is  computed  as  the  total  length 
of  the  segment  divided  by  the  number  of  points  to  be  plotted.  The  coordinates 
of  the  initial  node  are  stored  in  variables  XX  and  YY.  The  number  of  the  strip 
of  the  map  into  which  the  node  falls  is  computed.  Both  a current  value  of  the 
strip  number,  NMAP,  and  a value  for  the  previous  point,  NMAP0,  will  be  used. 

The  pen  position,  up  or  down,  is  determined  by  whether  the  initial  point  is  in 
bounds.  Variable  IPEN  will  be  3 if  the  point  is  out  of  bounds  and  2 if  the 
point  is  in  bounds.  If  the  point  is  out  of  bounds,  control  transfers  to  state- 
ment 130.  If  not,  the  coordinates  of  tl.a  point  are  converted  to  plotter  inches 
and  stored  in  variables  XP  and  YP. 
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Statement  130  starts  a loop  through  statement  170  that  will  advance 
the  pen  through  the  remaining  points  on  the  segment.  The  cumulative  length  is 
incremented  by  OS.  Subroutine  COORD  is  called  to  obtain  the  coordinates  of  the 
point  in  MCU. 

At  statement  140  the  coordinates  are  converted  to  plotter  inches.  The 
point  is  assumed  to  be  in  bounds,  and  variable  INB  is  set  to  1 . If  the  coord- 
inates of  the  point  are  out  of  bounds,  INB  is  reset  to  0.  If  the  pen  has  been 
up  and  the  current  point  is  out  of  bounds,  or  if  the  strip  number  is  greater 
than  the  number  of  the  final  strip,  control  transfers  to  statement  160.  Other- 
wise, the  pen  is  moved  to  the  position  of  the  current  point.  If  the  pen  is  up, 
it  is  lowered.  Variable  IPEN  is  recomputed  to  reflect  whether  the  point  is  in 
bounds. 


At  statement  160  the  number  of  the  current  strip  is  computed.  If  the 
current  strip  number  is  equal  to  the  previous  strip  number,  control  transfers 
to  statement  170.  If  not,  the  old  strip  number,  NMAP0,  is  set  equal  to  the 
current  strip  number.  IPEN  is  set  to  3 to  indicate  that  the  pen  is  up.  Con- 
trol transfers  to  statement  140,  where  the  coordinates  of  the  current  point  on 
the  new  strip  will  be  computed. 

Statement  170  is  the  end  of  the  loop  that  causes  the  segment  to  be 
drawn.  Statement  200  is  the  end  of  the  loop  that  draws  one  side  of  the  various 
segments.  Statement  210  marks  the  end  of  the  loop  on  each  side  of  the  street. 

Subroutine  SYMBOL  is  called  to  append  the  problem  title  to  the  lower 
left  corner  of  the  map.  A solid  line  is  drawn  2 inches  back  from  the  lower 
right  corner  by  two  calls  to  PLOT.  Subroutine  SYMBOL  appends  the  legend  TRAVEL 
to  the  line. 

The  loop  through  statement  230  controls  the  plotting  of  two  more  lines 
and  legends.  Subroutine  SYMBOL  plots  the  legend  COLLECT  BOTH  SIDES  on  the  first 
pass  through  the  loop  and  the  legend  COLLECT  RIGHT  SIDE  on  the  second  pass.  The 
loop  through  statement  220  plots  a horizontal  dashed  line.  On  the  second  pass 
through  the  230  loop,  small  vertical  lines  are  also  plotted. 


At  statement  300  the  plotter  pen  is  positioned  2 inches  beyond  the 
end  of  the  last  map  strip.  Control  returns  to  the  calling  program. 

m.  Subroutine  PATHPLT 

Subroutine  PATHPLT  draws  the  vehicle  path,  using  a solid  line  to 
indicate  travel  and  a dashed  line  to  indicate  collection.  The  subroutine  has 
six  arguments.  The  first  two  arguments,  NSP  and  NNP,  are  arrays  giving  the 
segment  and  node  numbers  in  the  path.  ,ne  third  argument,  CORT,  is  an  array 
of  collection-or-travel  indicators.  The  fourth,  NSG,  is  an  array  giving 
counts  of  segments  in  each  of  the  four  pieces  of  the  trip.  The  fifth  argu- 
ment, NTRIP,  is  either  1 or  2 depending  on  whether  the  trip  started  at  the 
garage  or  at  the  landfill.  The  sixth  argument,  NTPS,  gives  the  maximum  num- 
ber of  trips  per  day.  The  first  three  arguments  are  double-subscripted 
arrays.  The  first  subscript  corresponds  to  a step  in  the  path;  the  second 
corresponds  to  a piece  of  the  trip.  The  pieces  are  the  path  from  the  garage 
or  landfill  to  the  section;  the  path  within  the  section;  the  path  from  the 
section  to  the  landfill;  and,  for  the  final  trip  of  the  day,  the  path  from 
the  landfill  back  to  the  garage. 

A counter,  HAST,  is  initially  set  to  0.  The  count  of  trip  pieces, 
JF,  is  set  to  4.  If  the  current  trip  is  not  the  last  trip  of  the  day,  JF  is 
reset  to  3.  JF  is  assigned  the  value  3 in  the  next  statement.  This  state- 
ment causes  the  current  version  of  the  program  to  suppress  the  plotting  of 
the  trip  from  the  landfill  to  the  garage  and  can  be  deleted  if  the  user  wishes 
to  show  this  part  of  the  route. 

The  loop  through  statement  20  scans  each  piece  of  the  trip.  The 
number  of  segments  in  the  piece  is  saved  in  variable  N.  The  loop  on  state- 
ment 10  saves  the  segment  numbers  in  arrays  TRV  and  ISEG.  The  TRV  array  is 
sorted  into  increasing  order  by  subroutine  SHLSRT.  The  ISEG  array  is  carried 
along. 


The  loop  through  statement  40  removes  duplicate  segment  numbers 
from  the  ISEG  array  and  closes  up  the  empty  spaces.  Array  ITRV  is  set  to  0 
in  this  loop. 
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The  map  bounds  are  retrieved  from  variables  in  COMMON  block  MAPDATA 
and  are  stored  in  variables  XL,  XR,  YB,  and  YT.  The  maximum  length,  XMX;  the 
number  of  map  strips,  MX;  the  map  scale  factors;  and  the  intervals  in  MCU  at 
which  the  strips  are  to  be  cut  are  computed. 
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The  outer  loop  through  statement  160  controls  the  scanning  of  the 
pieces  of  the  trip.  The  inner  loop  through  statement  160  will  plot  each  path 
segment  that  falls  within  the  frame  of  the  map.  Variables  NI  and  NF  are  set 
equal  to  the  numbers  of  the  nodes  bounding  the  path  segment.  The  midpoint  co- 
ordinates of  the  segment  are  saved  in  variables  XMD  and  YMD.  The  lines  in  the 
node-number  array  at  which  the  nodes  occur  are  saved  in  variables  LI  and  LF. 

If  the  nodes  have  been  found  in  the  NODNUM  array,  control  transfers  to  state- 
ment 60.  Otherwise,  an  error  message  is  printed  and  control  transfers  to  state- 
ment 160.  At  statement  60,  the  node  numbers  are  compared  with  the  numbers  of 
the  nodes  bounding  the  segment  in  the  segment  data.  If  the  numbers  are  equal, 
control  transfers  to  statement  80.  Otherwise,  an  error  message  is  printed  and 
control  transfers  to  statement  160. 

Following  statement  80,  the  coordinates  of  the  nodes  are  retrieved. 

The  segment  is  assumed  to  be  entirely  within  the  bounds  of  the  map,  and  indi- 
cators INBI,  INBM,  and  INBF  are  set  to  1 . If  the  coordinates  of  the  initial 
node  lie  outside  the  frame  of  the  map,  INBI  is  set  to  0.  Similar  tests  are 
made  on  the  coordinates  of  the  midpoint  and  final  nodes  of  the  segment.  If 
all  three  points  are  outside  the  frame  of  the  map,  control  transfers  to  state- 
ment 160  and  no  path  is  plotted  for  that  segment.  For  segments  that  are  at 
least  partially  within  the  frame,  the  total  length  of  the  segment,  in  miles, 
is  stored  in  variable  T0TLEN.  If  the  segment  is  straight,  T0TLEN  is  computed 
using  the  end-point  coordinates  and  the  map  distance  conversion  factor.  The 
number  of  1/10-inch  steps  to  the  middle  of  the  segment,  NPMID,  is  computed. 

If  NPMID  is  0,  it  is  reset  to  1.  The  number  of  points  per  segment,  NPPSEG,  is 
one  less  than  twice  NPMID.  Variable  CUMLEN  is  set  to  0.  A step  size,  DS,  is 
computed  from  the  total  length  and  the  number  of  points  per  segment. 

Variable  ISH  is  set  to  either  1 or  8,  depending  on  the  direction  of 
travel  on  the  segment.  Similarly,  variable  DIR  is  set  to  1 or  -1.  Segment  KK 
is  found  in  the  ISEG  array,  and  the  line  number  is  stored  in  variable  LTR.  The 
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number  of  times  the  segment  has  been  traversed,  NTRV,  is  retrieved  from  the 
ITRV  array.  In  the  ITRV  array,  the  octal  units  digit  gives  the  number  of  tra- 
versals from  starting  to  ending  node,  and  the  octal  tens  digit  gives  the  number 
of  traversals  in  the  other  direction.  The  distance  of  the  path  line  from  the 
center  of  the  street,  W,  is  computed  on  the  basis  of  the  number  of  previous 
traversals.  The  appropriate  digit  in  the  ITRV  array  is  incremented  by  ISH. 

The  col lection-or- travel  indicator  is  retrieved  and  stored  in  variable  ACT. 
Logical  variable  RSO  is  computed  to  indicate  whether  collection  is  from  the 
right  side  only. 

Subroutine  SHAPCOM  is  called  to  set  up  the  parameters  needed  to  gen- 
erate coordinates  of  the  points  on  the  path.  Subroutine  COORD  is  called  to  ob- 
tain the  coordinates  of  the  starting  point  of  the  segment,  in  MCU.  The  number 
of  the  strip  of  the  map,  NMAP,  is  computed.  The  pen  position,  up  or  down,  is 
determined  by  whether  the  initial  point  was  in  bounds.  The  computation  for 
variable  IPEN  yields  3 if  the  point  is  out  of  bounds  or  2 if  the  point  is  in 
bounds.  If  the  point  is  out  of  bounds  or  if  an  error  has  been  detected  by  sub- 
routine COORD,  control  transfers  to  statement  100.  Otherwise,  the  coordinates 
of  the  point  are  saved  in  variables  XLAST  and  YLAST.  The  coordinates  are  con- 
verted to  plotter  inches  and  stored  in  variables  XP  and  YP.  If  the  current 
point  is  the  first  to  be  plotted  or  if  the  current  strip  number  differs  from 
the  previous  number,  subroutine  PLOT  is  called  to  raise  the  pen  and  move  it  to 
the  current  point.  Variable  IFIRST  is  set  to  0.  Variable  NMAPO  is  set  equal 
to  the  current  strip  number.  The  pen  is  moved  to  the  current  point  in  the  down 
position. 


Statement  100  starts  a loop  through  statement  150  that  will  advance 
the  pen  through  the  remaining  points  on  the  segment.  The  cumulative  length  is 
incremented  by  DS.  A count  of  points  drawn,  KT0T,  is  incremented  by  1.  Sub- 
routine COORD  is  called  to  obtain  the  coordinates  of  the  point  in  MCU.  If 
COORD  finds  an  error,  control  transfers  to  statement  150. 

At  statement  110  the  coordinates  of  the  previous  point  are  saved  in 
variables  XLAST  and  YLAST.  The  coordinates  of  the  current  point  are  computed 
in  plotter  inches  and  stored  in  variables  XP  and  YP.  The  point  is  assumed  to 
be  in  bounds,  and  variable  INB  is  set  to  1 . If  the  coordinates  of  the  point 
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are  out  of  bounds,  INB  is  reset  to  0.  If  the  path  is  to  indicate  collection, 
control  transfers  to  statement  120.  Otherwise,  if  the  pen  is  in  the  up  posi- 
tion, it  is  lowered.  Variable  IPEN  is  recomputed  to  reflect  whether  the  point 
is  in  bounds.  If  the  count  of  points,  KTOT,  is  a multiple  of  KARO  (in  the  cur- 
rent version  of  the  program  KARO  is  equal  to  30),  control  transfers  to  state- 
ment 125.  If  not,  control  transfers  to  statement  140. 

At  statement  120,  IPEN  is  set  to  2.  If  the  loop  index,  K,  is  even, 
control  transfers  to  statement  130.  Otherwise,  IPEN  is  set  equal  to  3.  If 
KTOT  modulo  KARO  is  greater  than  1,  control  transfers  to  statement  140.  Other- 
wise, starting  at  statement  125,  displacements  are  computed  for  drawing  an 
arrowhead.  The  arrowhead  is  drawn  by  three  calls  to  PLOT.  Control  transfers 
to  statement  140. 

At  statement  130,  if  collection  is  not  from  the  right  side  only,  con- 
trol transfers  to  statement  140.  Otherwise,  a small  mark  is  drawn  perpendicular 
to  the  direction  of  the  path.  At  statement  140,  the  strip  number  of  the  point 
is  computed.  If  the  strip  number  equals  the  previous  strip  number,  control 
transfers  to  statement  150.  Otherwise,  if  KTOT  modulo  KARO  is  less  than  or 
equal  to  1,  KTOT  is  incremented  by  2.  NMAPO  is  set  equal  to  the  current  strip 
number.  IPEN  is  set  equal  to  3.  Control  returns  to  statement  110,  where  the 
pen  will  be  positioned  at  the  current  point  on  the  new  strip. 

Statement  150  marks  the  end  of  the  loop  that  plots  the  points  on  the 
current  path  segment.  Statement  160  marks  the  end  of  the  two  loops  on  the  seg- 
ments in  the  pieces  of  the  travel  path.  Control  returns  to  the  calling  program. 

n.  Program  PHASE4 

Program  PHASE4  uses  the  data  prepared  by  programs  RCINPT  and  PHASE3 
to  select  a refuse-collection  schedule  for  a particular  area,  prints  that  sched- 
ule, and  plots  a map  showing  the  collection  routes.  The  file  assignments  in- 
clude TAPE1,  a binary  file  containing  segment  data;  TAPE2,  another  binary  file 
containing  node  data;  TAPE3,  street-name  data  buffered  in  for  use;  TAPE8,  the 
Calcomp  plot  tape;  TAPE9,  the  formatted  path  data;  and  TAPE5,  the  input  file. 
COMMON  blocks  are  set  up  and  variables  dimensioned.  A title,  the  unit  for  the 


measurement  of  refuse  quantity,  and  an  indicator  that  specifies  whether  the 
trips  are  to  be  paired  are  read  from  cards. 


The  segment  data  are  then  read  from  TAPE1 . These  data  include  the 
number  of  segments  and  the  following  specific  information  for  each  segment: 

NSTR,  the  street  number;  NN1,  the  initial  node  for  the  segment;  NN2,  the  final 
node  for  the  segment;  FLEN,  the  length  of  the  segment;  NH,  the  number  of  houses 
on  the  segment;  FMPH,  the  speed  limit  on  the  segment,  in  miles  per  hour;  NWAY, 
which  indicates  whether  the  vehicle  will  be  on  a one-way  or  a two-way  street; 
RQF,  the  refuse-quantity  adjustment  factor;  XMID,  the  x-coordinate  of  the  mid- 
point of  the  segment;  YMID,  the  y-coordinate  of  the  midpoint  of  the  segment; 
and  SF,  the  shape  code.  AVMD,  the  map  distance  conversion  factor,  is  the  last 
word  on  the  file. 

The  node  data  are  then  read  from  file  TAPE2.  The  data  include  NHTOT, 
the  total  number  of  houses;  TOTREF,  the  total  amount  of  refuse;  KNODES,  the 
total  number  of  nodes;  and  the  following  specific  data  for  each  node:  NODNUM, 
the  number  of  the  node;  NBS,  the  numbers  of  segments  bounding  the  node;  XNOD, 
the  x-coordinate  of  the  node;  and  YNOD,  the  y-coordinate  of  the  node. 

The  title  is  then  printed,  along  with  the  count  of  street  segments 
and  the  count  of  nodes.  Subroutine  STRIN  is  called  to  read  and  save  the  street- 
name  information.  The  following  time-limitation  data  are  read  from  cards: 
TSTOPH,  the  time  spent  at  each  house;  TSTOPR,  the  time  required  to  collect  each 
unit  of  refuse;  TUNLD,  the  time  required  to  unload  the  vehicle;  TMXTR,  the  maxi- 
mum time  allowed  for  each  vehicle  trip;  TMXDAY,  the  maximum  time  allowed  for 
each  vehicle  to  be  in  operation  during  a day;  TSTART,  the  starting  time  of  the 
work  day;  SCOLL,  the  speed  of  the  vehicle  during  refuse  collection;  DLIINCH,  the 
duration  of  the  lunch  break;  TLUNCH,  the  time  the  lunch  break  should  begin; 
array  DBRK,  the  durations  of  the  mid-morning  and  mid-afternoon  breaks;  and  array 
TBRK,  the  times  at  which  the  mid-morning  and  mid-afternoon  breaks  should  start. 
Some  of  the  data  are  then  printed  out  as  a check  on  the  input. 


The  maximum  trip  time,  the  maximum  working  time  per  day,  the  starting 
time  for  the  day,  and  the  vehicle  speed  during  collection  are  then  examined. 


Each  value  read  should  be  greater  than  zero.  If  the  value  of  any  of  these  vari- 
ables is  not  greater  than  zero,  PHASE4  will  change  it  to  a predetermined  value. 
TMXTR  will  be  set  to  4 hours,  TMXDAY  to  8 hours,  TSTART  to  8.00  hours,  and  SCOLL 
to  5 miles  per  hour.  The  rest  of  the  input  data  are  then  printed. 

The  vehicle  capacities  and  identifications  are  read  next.  TC,  the 
truck  capacity , and  VID,  the  vehicle  identification,  are  read  for  up  to  ten 
truck  capacities  with  five  words  per  vehicle  identification.  After  they  have 
been  read,  the  data  are  printed.  The  loop  is  executed  up  to  11  times,  as 
necessary,  to  reach  the  end-of-record  card  terminating  the  data. 

Program  PHASE4  next  checks  for  map-bounds  cards.  The  data  read  from 
cards  are  NSCN,  the  section  number;  NTRP,  the  number  of  the  trip;  XMN,  the  min- 
imum x-value  for  the  map;  XMX,  the  maximum  x-value;  XLN,  the  length  in  the  x- 
direction;  YMN,  the  minimum  y-value;  YMX,  the  maximum  y-value;  and  YLN,  the 
length  in  the  y-direction.  If  no  map  bounds  are  specified  for  a particular 
map,  the  map  will  show  travel  in  the  collection  region  but  not  necessarily  the 
path  to  or  from  the  garage  or  landfill.  If  bounds  are  specified  for  any  of  the 
maps,  the  data  are  printed  out. 

The  path  data  are  read  from  TAPE9  for  the  travel  from  the  landfill 
back  to  the  garage.  These  data  include  N,  the  number  of  segments  in  the  path; 
DIST,  the  distance  covered;  NNP,  the  node  numbers;  NSP,  the  segment  numbers; 
CORT,  the  collection-or-travel  indicators;  and  the  final  node  number.  If  the 
path  data  are  not  found  on  TAPE9,  PHASE4  will  check  for  data  on  cards  on  input 
file  TAPE5.  If  the  data  are  on  TAPE5,  TAPE9  is  rewound  and  the  card  images  are 
written  to  file  TAPE9.  If  no  path  data  are  found  on  TAPE9  or  on  cards,  the  job 
is  terminated. 

After  the  data  files  have  been  examined,  plotting  is  initialized  by 
calls  to  PLOTS  and  PLOT.  The  first  call  to  CUMTD  determines  the  time  and  dis- 
tance from  the  landfill  to  the  garage.  The  following  variables  are  now  set  to 
0:  NSCOLD,  the  old  section  number;  NTRIP,  the  trip  number;  TOTD,  the  total  dis- 
tance; TOTT,  the  total  time;  and  TOTR,  the  total  refuse  collected.  Data  for 
travel  to,  within,  and  from  the  collection  region  are  read.  If  the  data  are 
entered  from  cards  (TAPE5),  they  are  written  to  file  TAPE9.  The  data  consist 
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of  N,  the  number  of  segments  in  the  path;  DIST,  the  distance  covered;  NSC,  the 
section  number;  TRC,  the  truck  capacity;  TRL , the  refuse  load;  NNP,  the  numbers 
of  the  nodes  in  the  path;  NSP,  the  numbers  of  the  segments  in  the  path;  and 
CORT,  the  collection-or-travel  indicators.  A call  to  CUMTD  provides  the  total 
time,  distance  covered,  and  refuse  collected  for  each  part  of  the  trip.  These 
totals  are  accumulated  in  variables  TOTD  (total  distance),  TOTT  (total  time), 
and  TOTR  (total  refuse  collected).  If  the  part  of  the  trip  just  completed  was 
travel  within  the  collection  area,  the  number  of  houses  serviced  is  stored  in 
variable  NHS. 

PHASE4  next  determines  whether  map-bounds  cards  have  been  read  for 
the  collection  region  just  processed.  If  so,  the  program  continues  at  state- 
ment 480.  If  not,  XMAX,  XMIN,  YMAX,  and  YMIN  are  set  to  default  values.  These 
default  values  are  determined  by  a consideration  of  both  the  nodes  ar.d  the  seg- 
ments within  the  collection  region.  The  section  number  and  the  initial  values 
of  XMAX,  XMIN,  YMAX,  and  YMIN  are  printed  out.  A scaling  factor,  SC,  is  deter- 
mined, and  the  values  are  adjusted  in  preparation  for  plotting  the  map  of  the 
collection  region.  The  adjusted  values  are  then  printed. 

If  path  data  are  read  from  cards,  an  end-of-file  mark  is  placed  on 
TAPE9,  and  the  file  is  rewound.  The  trip  information,  consisting  of  the  sec- 
tion number,  trip  number,  distance  covered,  time,  number  of  households  serviced, 
capacity,  and  load,  is  printed. 

If  the  trips  are  not  to  be  paired  in  the  order  in  which  they  occur, 
subroutine  SHLSRT  is  called  to  sort  the  vehicles  by  capacity.  Pointers  to  the 
trips  in  the  IORD  array  are  carried  along  during  the  sort. 

In  the  loop  through  statement  540,  the  travel  times  for  each  vehicle 
capacity  are  sorted.  If  the  longest  time  exceeds  the  maximum  trip  time  speci- 
fied by  the  user,  a warning  message  is  printed  and  the  maximum  trip  time  is  ex- 
tended. The  maximum  trip  time  for  vehicles  of  each  capacity  is  saved  in  the 
TMXTRV  array. 

A page  heading  for  the  final  route  summary  is  printed.  The  loop 
through  statement  620  accumulates  the  values  printed  in  the  route  summary.  The 


loop  through  statement  550  seeks  the  line  number  of  the  vehicle  identification 
corresponding  to  the  vehicle  capacity  that  is  now  being  processed  by  loop  620. 

The  loop  through  statement  570  examines  the  trips  selected  for  vehi- 
cles of  the  capacity  now  being  processed.  The  distance  saved  by  using  the  sec- 
ond rather  than  the  first  trip  choice  for  each  section  is  stored  in  array  TEM. 

If  the  time  for  the  first  trip  choice  exceeds  the  maximum  trip  time,  the  dis- 
tance saving  is  increased  by  2000.  If  the  trip  time  for  the  second  choice  ex- 
ceeds the  maximum  trip  time,  the  distance  saving  is  decreased  by  2000.  This 
adjustment  is  made  so  that  when  one  trip  choice  exceeds  the  maximum  time  and 
the  other  does  not,  the  program  will  select  the  trip  that  takes  less  time. 

Subroutine  SHLSRT  is  called  to  sort  the  distance  savings  into  in- 
creasing order.  The  loop  on  statement  585  replaces  the  distance  savings  in 
array  TEM  by  the  time  for  the  appropriate  trips.  If  more  than  one  route  is 
present  and  if  there  are  two  choices  of  trip  per  section,  subroutine  SHLSRT 
is  called  to  sort  the  times  of  the  first  trips  into  increasing  order.  The 
trip  pointers  are  carried  along  during  the  sort.  If  there  is  more  than  one 
afternoon  trip,  subroutine  SHLSRT  is  called  to  sort  the  afternoon  trip  times 
into  increasing  order.  The  trip  pointers  are  carried  along  during  the  sort. 

The  loop  through  statement  610  accumulates  the  total  number  of  houses, 
refuse  quantity,  time,  and  distance  for  each  route.  These  values  are  printed 
as  part  of  the  route  summary.  The  section  numbers  of  the  trips  comprising  the 
route  are  saved  in  the  IRS  array.  Cumulative  values  of  distance,  time,  refuse, 
and  houses  serviced  are  also  accumulated  in  the  loop.  Following  statement  620, 
the  totals  for  all  routes  are  printed. 

File  TAPE9,  which  contains  the  route  data,  is  rewound.  The  sequence 
number  of  the  last  trip  read,  LSEQN,  is  set  to  -1  to  indicate  that  the  file  is 
positioned  at  the  beginning  of  the  data.  The  trip  from  the  landfill  to  the 
garage  has  sequence  number  0.  The  trips  servicing  each  section  have  sequence 
numbers  1 through  the  number  of  trips. 

The  loop  through  statement  680  generates  the  maps  and  schedule  for 
each  route.  The  loop  is  executed  one  time  more  than  the  number  of  routes  so 
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that  maps  without  paths  can  also  be  plotted.  Variable  JTRIPS  is  computed  to 
be  the  number  of  trips  per  route.  The  loop  through  statement  670  controls 
schedule  and  map  generation  for  each  trip.  The  section  number  is  obtained  from 
the  IRS  array.  The  sequence  number  for  the  trip,  ISEQN,  is  computed.  Subrou- 
tine P0SIT9  is  called  to  position  TAPE9  at  the  beginning  of  the  desired  trip. 
The  loop  through  statement  640  reads  the  three  pieces  of  the  trip  from  TAPE9. 

If  the  trip  is  the  first  trip  of  the  route,  the  route  number,  the 
problem  title,  and  the  vehicle  identification  are  printed  as  the  heading  for 
the  schedule.  Subroutine  PRSCHED  is  called  to  print  the  schedule  for  the  trip. 

The  loop  through  statement  660  examines  the  map-bounds  data.  The 
bounds  and  lengths  of  the  map  are  retrieved.  If  the  bounds  in  the  x or  y 
directions  allow  no  width  or  height  to  the  map,  the  map  plotting  is  cypassed 
and  control  transfers  to  statement  670.  If  the  loop  index  indicates  a route 
with  trips  in  it,  subroutine  PATHPLT  is  called  to  plot  the  path.  Subroutines 
NUMBER  and  SYMBOL  are  called  to  append  a legend  to  the  map.  Subroutine  STNAME 
is  called  to  append  street  names  to  the  map.  Subroutine  MAPPLT  is  called  to 
draw  the  sides  of  the  street  on  the  map. 

Statements  660,  670,  and  680  mark  the  ends  of  their  respective  loops. 

Subroutine  PLOT  is  called  to  terminate  the  plot  file.  Subroutine 
EXIT  is  called  to  return  control  to  the  system. 
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SECTION  IV 
INPUT  AND  OUTPUT 


1 . INPUT 

Input  to  program  PHASE4  consists  of  card  input  and  four  disk  files.  Three 
of  the  disk  files  are  generated  by  program  RCINPT  and  one  by  program  PHASE3. 

The  files  are  read  by  the  main  program,  PHASE4,  and  by  subroutine  STRIN. 

a.  Card  Input 

The  form  and  contents  of  the  data  cards  are  shown  in  Table  1.  Four 
types  of  data  cards  are  required  for  program  PHASE4;  three  more  types  are  op- 
tional. The  required  cards  are  a title  card,  a refuse-unit  card,  two  time- 
restriction  cards,  and  vehicle-identification  cards.  The  optional  cards  are 
map-bounds,  street-name,  and  path-information  cards. 

Data  cards  for  Kirtland  Air  Force  Base  are  shown  in  Appendix  D.  A 
record  of  map-bounds  cards  is  included.  In  this  case,  there  are  no  street-name 
or  path-information  cards  because  these  data  are  obtained  from  disk  files.  A 
detailed  description  of  the  preparation  of  these  data  cards  can  be  found  in 
Reference  1 . 

b.  Disk  Files 

Disk  file  TAPE1  contains  segment  data  and  the  map  distance  conversion 
factor  (miles  per  MCU)  for  the  overall  map.  All  of  the  data  are  read  by  one 
binary  READ  statement.  The  first  word  is  the  count  of  the  segments.  The  seg- 
ment data  follow,  11  words  per  segment.  After  the  segment  data  comes  the  over- 
all distance  conversion  factor.  The  following  list  is  used  in  the  READ  state- 
ment: 

Reference 

'luzzolino,  Harold  J.,  Air  Force  Refuse-Collection  Scheduling  Program , CEEDO-TR 
77-32,  Tyndall  Air  Force  Base,  Florida,  January  1978. 
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TABLE  1.  PHASE4  DATA  CARDS 


Record  Card  Columns  Format 


Contents 


2 1-20  2A10  Refuse  units 

21-25  15  Trip  pairing  control 

The  following  card  is  optional  and  may  be  repeated  up  to  300  times. 

3 1-5  15  Street  number  (right  justified) 

11-60  5A10  Street  name  (left  justified) 

End  of  Record  7-8-9  multipunched  in  column  1. 


If  card  3 is  omitted,  omit  End  of  Record  and  append  the  next  three 
cards  to  record  1 . I 


11-20 

21-30 


F10.0 

F10.0 

F10.0 


31-40  F10.0 


41-50 

51-60 


F10.0 

F10.0 


61-70  F10.0 


1-10  F10.0 

11-20  F10.0 


21-30 

31-40 

41-50 


F10.0 

F10.0 

F10.0 


51-60  F10.0 


Stop  time  per  household,  in  minutes 

Stop  time  per  unit  refuse,  in  minutes 

Unloading  time,  in  minutes 

Maximum  trip  time,  in  hours 

Maximum  working  time  per  day,  in  hours 

Starting  time,  in  hours 

Average  vehicle  speed,  in  mph,  between  two 
collection  points  on  a street 

Duration  of  lunch,  in  minutes 
Starting  time  of  lunch,  in  hours 
Duration  of  first  break,  in  minutes 
Starting  time  of  first  break,  in  hours 
Duration  of  second  break,  in  minutes 
Starting  time  of  second  break,  in  hours 


The  following  card  may  be  repeated  up  to  10  times. 


« 

End  of  Record 


1-10  F10.2  Vehicle  capacity 

11-60  5A1C  Vehicle  identification 

7-8-9  multipunched  in  column  1. 
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TABLE  1.  PHASE4  DATA  CARDS  (Concluded) 


Record  Card  Columns  Format 


Contents 


The  following  card  is  optional  and  may  be  repeated  up  to 
100  times.  I I 


1-5 

15 

6-10 

15 

11-20 

F10.0 

21-30 

F10. 0 

31-40 

F10.0 

41-50 

F10.0 

51-60 

F10. 0 

61-70 

F10.0 

End  of  Record 


Section  number 
AM/PM  trip  indicator 
Minimum  x-coordinate 


Minimum  y-coordinate 
Maximum  y-coordinate 
Length  of  y-di recti  on,  in  inch< 
7-8-9  multipunched  in  column  1 


The  following  record  is  used  only  if  path  data  are  not 
on  file  TAPE9.  I 1 


Count  of  segments  in  this  piece 

Total  distance  for  this  piece  of  the  trip, 
in  miles 

Section  number 

Vehicle  capacity 

Total  refuse  quantity  for  this  trip 

Node  number 
Segment  number 

Col lection-or-travel  indicator 


1-5 

15 

6-10 

F10.3 

11-15 

15 

16-25 

F10.3 

26-35 

F10.3 

1-5 

15 

6-9 

14 

10 

A1 

The  above  three  items  may  be  repeated  to  the  end  of  the  card. 
This  card  is  repeated  until  the  entire  path  to  the  collection 
region  has  been  described., 


The  two  cards  are  used  twice  more  for  each  trip,  once  to  de- 
scribe the  path  in  the  collection  region  and  again  to  describe 
the  path  from  the  collection  region  to  the  landfill. 


7-8-9  multipunched  in  column  1. 


End  of  Record 
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1 
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NSEG,  (NSTR(I),  NN1(I),  NN2(I),  FLEN(I),  NH(I),  FMPH(I), 
NWAY(I),  RQF(I),  XMID(I),  YMID(I),  SF(I),  I = 1,  NSEG),  AVMD 

The  items  in  the  list  are  defined  as  follows: 


NSEG 

NSTR(I) 

NN1(I) 

NN2( I ) 

FLEN(I) 

NH  ( I ) 

FMPH(I) 

NWAY(I) 

RQF(I) 

XMID(I) 

YMID(I) 

SF  ( I ) 

AVMD 


number  of  segments  

street  number 
starting  node  number 
ending  node  number 
street  length,  in  miles 

number  of  houses  

speed  limit,  in  mph 
number  of  ways  of  travel 
refuse-quantity  adjustment  factor 
x-coordinate  of  segment  midpoint 
y-coordinate  of  segment  midpoint 
shape  code 

map  distance  conversion  factor,  in  miles  per  MCU 


- These  are  repeated 
for  each  segment. 


The  sign  of  the  number  of  houses  is  used  to  indicate  whether  collection  is  on 
both  sides  or  on  only  the  right  side  of  the  segment.  NH  is  negative  to  indi- 
cate collection  on  only  the  right  side  of  the  street,  or  positive  to  indicate 
collection  from  both  sides  of  the  street. 


Disk  file  TAPE2  contains  refuse-quantity  information  and  the  node 
data.  All  of  the  data  are  read  by  one  binary  READ  statement.  The  first  three 
words  are  the  total  number  of  houses  or  stops,  the  total  refuse  quantity,  and 
a count  of  the  nodes.  The  node  data  follow,  four  words  per  node.  The  follow- 
ing list  is  used  in  the  READ  statement: 

NHTOT,  TOTREF,  KNODES,  (NODNUM(I),  NBS(I),  XNOD(I),  YNOD(I),  1=1,  KNODES) 

The  items  in  the  list  are  defined  as  follows: 

NHTOT  = total  number  of  houses  or  stops 
TOTREF  = total  refuse  quantity 

KNODES  = count  of  nodes 


» 
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NODNUM(I)  = node  number 

NBS(I)  = packed  bounding-segment  numbers 

XNOD(I)  = x-coordinate  of  node 

YNOD(I)  = y-coordinate  of  node 


The  index  I corresponds  to  the  Ith  node.  Variable  NBS(I)  contains  up  to  6 seg- 
ment numbers,  each  occupying  10  of  the  60  bits,  for  segments  bounding  node 
NODNUM(I).  Since  node  numbers  are  assigned  by  the  user,  NODNUM(I)  usually  is 
not  the  same  as  I. 


- These  are  repeated 
for  each  node. 


Disk  file  TAPES  contains  the  street  numbers  and  names  (70  characters 
each).  The  data  are  read  100  streets  at  a time,  using  a BUFFER  IN  statement. 

Up  to  300  streets  may  be  read.  The  data  consist  of  all  of  array  NUMSTR  (street 
numbers),  followed  by  all  of  array  NAMSTR  (street  names).  Because  the  arrays 
are  adjacent  in  storage,  each  record  appears  as  follows: 

NUMSTR(l)  = number  of  first  street 


NUMSTR(IOO) 
NAMSTR (1 ,1 ) 

NAMSTR (7,1 ) 

NAMSTR (1 ,2)  = 


number  of  100th  street 
7-word  name  of  first  street 

first  word  of  second  street  name 


NAMSTR (7, 100) 


last  word  of  100th  street  name 


The  last  record  contains  zeros  in  unused  words. 


Disk  file  TAPE9  contains  the  path  data.  A description  of  the  path 
from  the  landfill  to  the  garage  occurs  at  the  beginning  of  the  file.  Two  trips 
follow  for  each  section,  one  starting  at  the  garage  and  one  starting  at  the 
landfill.  Each  trip  is  divided  into  three  pieces,  and  a header  precedes  each 
piece. 


The  path  from  the  landfill  to  the  garage  and  its  header  are  read  at 
statement  270  in  the  main  program  using  a formatted  READ  statement.  The  re- 
maining trips  and  their  headers  are  read  in  the  loop  through  statement  420  in 
the  main  program.  The  following  list  is  used  in  reading  the  trip  pieces: 
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N, DIST( J), NSC, TRC,TRL,(NNP(K,J),NSP(K,J),C0RT(K,J),K=1 ,N)  ,NNP(N+1 ,J) 


The  subscript  J has  values  1 , 2,  or  3 for  the  pieces  of  the  trips  servicing 
sections  and  value  4 for  the  trip  from  the  landfill  to  the  garage.  The  items 
in  the  list  are 


N 

DIST(J) 

NSC 

TRC 

TRL 

NNP(K, J) 
NSP(K.J) 
C0RT(K,J) 
NNP(N+1,J) 


Count  of  segments  in  piece 
Length  of  piece,  in  miles 
Section  number 
Vehicle  capacity 

Vehicle  load  

Node  number 
Segment  number 

Col lection-or- travel  indicator 


These  data  comprise  the  header. 


Jhese  data  are  repeated  for 
‘each  segment  in  the  path. 


Number  of  the  node  terminating  the  trip 


The  trip  from  the  landfill  to  the  garage  lacks  items  NSC,  TRC,  and 
TRL  in  its  header.  The  format  for  its  header  is  I5.F10.0.  The  format  used 
for  the  remaining  headers  and  path  data  is  I5,F10.3,I5,2F10.3/8( 15,14 ,A1 ) . 

The  route  data  are  read  from  either  TAPE9  or  cards.  If  the  data  are  on  cards, 
they  are  written  to  TAPE9  in  the  above  format. 

TAPE9  is  reread  in  the  loop  through  statement  640  in  the  main  pro- 
gram. At  this  time  the  data  are  used  to  generate  the  schedules  and  maps. 


2.  OUTPUT 

a.  Plot  File 

File  TAPE8,  the  plot  file,  will  be  on  disk  or  tape  depending  on  the 
procedure  used  by  the  local  installation  to  produce  plots.  Each  map  occupies 
one  file. 


The  final  maps  show  travel  to  the  collection  region,  within  the  col- 
lection region,  and  from  the  collection  region  to  the  landfill.  Street  names 
are  appended  to  the  map,  and  the  direction  of  travel  is  shown  by  arrowheads  on 
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the  lines  indicating  collection  or  travel.  Subsequent  traversals  along  the 
same  street  segment  are  indicated  by  path  lines  drawn  closer  to  the  center  of 
the  street.  Maps  for  Kirtland  Air  Force  Base  are  given  in  Appendix  E.  The 
maps  were  produced  using  the  map-bounds  cards  in  Appendix  0 and  have  been  re- 
duced to  half  size. 

b.  Printed  Output 

The  printed  output  consists  of  five  sections:  a listing  of  input 
data,  map-bounds  listings,  a trip  summary,  a final  route  summary,  and  detailed 
descriptions  of  the  routes. 

The  listing  of  input  data  is  provided  for  verification  of  the  card 
input.  The  problem  title  is  printed  first.  Counts  of  the  segments,  nodes,  and 
cards  containing  street-names  are  printed.  Also  listed  are  the  unit  of  refuse 
and  the  time  restrictions. 

If  bounds  were  specified  for  any  of  the  maps,  these  data  are  printed 
on  the  second  page  of  output.  The  third  page  lists  map  bounds  selected  by  the 
program  for  sections  not  given  map  bounds  by  the  user.  The  column  headings 
are  printed  even  if  the  program  selects  no  bounds. 

The  trip  summary  gives  general  information  about  the  trip  choices  for 
each  section. 

A final  route  summary  is  printed  on  the  fifth  page  of  output.  The 
vehicle  identification,  vehicle  capacity,  section  numbers  of  the  trips,  dis- 
tance, time,  households  serviced,  and  refuse  quantity  are  listed  for  each 
route.  Totals  are  given  for  the  distance,  time,  households  serviced,  and 
refuse  quantity  for  all  the  routes. 

The  detailed  route  description  shows  the  action  to  be  taken  on  each 
street  of  each  route.  The  first  action  printed  instructs  the  driver  to  leave 
the  garage.  Then  instructions  are  given  either  to  drive  on  or  collect  from 
the  street  up  to  a particular  cross  street.  The  speed  limit,  distance  traveled, 
and  number  of  households  serviced  on  that  street  are  given.  An  estimate  of  the 


41 


r 


„ I,  collection  is 

t,0„  is  to  be  comP'e“d5  print!  *«er  ““T 

. « * *h'clTotaU^  e u^\  «• 

fled,  the  Perce  reqion,  the  p*  A unloading  o 

* “wt^  *«»  * *•  *7 

is— . rW€iW.^— 


U F- 


time  of  day  at  which  the  action  is  to  be  completed  is  given.  If  collection  is 
performed,  the  percent  of  vehicle  capacity  used  is  printed.  After  the  descrip- 
tion of  collection  in  the  region,  the  path  to  the  landfill  is  described.  At 
the  end  of  the  description  of  the  day's  final  unloading  of  the  vehicle,  the 
path  back  to  the  garage  is  described. 

Complete  printed  output  for  Kirtland  Air  Force  Base  is  given  in 
Appendix  F. 


1 
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SECTION  V 

PROGRAM  REQUIREMENTS 


) 


1 . SYSTEM 

Program  PHASE4  is  written  entirely  in  FORTRAN  IV.  The  program  runs  on  a 
CDC  6600  computer  using  a SCOPE  3.4.4  operating  system. 

Ten  obvious  types  of  computer-dependent  coding  occur  in  program  PHASE4  and 
its  subroutines.  Subroutine  PRSCHED  assumes  a 130-character  output  line.  A 
60-bit  word  is  assumed  in  subroutine  NUMBER.  Subroutine  PRSCHED  assumes  ten 
characters  per  word.  System  subroutine  SHIFT  is  used  by  subroutines  NUMBER  and 
PRSCHED.  An  ENCODE  statement  is  used  in  subroutine  NUMBER  and  function  HM;  a 
DECODE  statement  is  used  in  subroutine  STNAME . The  AND  masking  operation  is 
used  in  subroutines  NUMBER,  PRSCHED,  and  PATHPLT.  Asterisk-bounded  text  is 
used  in  format  statements  in  the  main  program  and  subroutines  PATHPLT,  PRSCHED, 
and  STRIN.  Multiple  replacement  statements  occur  in  the  main  program  and  in 
subroutines  MAPPLT,  SHAPCOM,  CUMTD,  and  NUMBER.  A dollar  sign  is  used  to  sep- 
arate FORTRAN  statements  in  the  main  program  and  in  all  subprograms  except 
function  I FIND  and  subroutines  STRIN  and  NUMBER.  An  R-format  text  variable  is 
used  in  subroutines  SHAPCOM  and  COORD. 

More  subtle  types  of  machine  dependencies  may  exist,  according  to  the 
machine  used. 


2.  STORAGE 

The  core  requirement  is  slightly  less  than  111  ,0008  words.  If  data  are 
written  to  file  TAPE9,  the  maximum  peripheral  storage  used  by  the  file  will  not 
exceed  33,800  words.  The  maximum  storage  required  by  the  plot  file,  TAPE8, 
should  not  exceed  30,000,000  words,  although  an  estimate  of  50,000  words  per 
map  is  more  typical. 


43 


3.  TIME 


The  running  time  for  PHASE4  varies  with  the  number  of  maps  and  the  number 
of  segments  per  map.  Since  the  number  of  segments  per  map  can  vary  consider- 
ably, the  following  figures  are  estimates: 

CP  time  = 2 seconds  +2.5  seconds  x (number  of  maps) 

I/O  time  = 2 seconds  +2.5  seconds  x (number  of  maps) 

PP  time  = 30  seconds  + 3 seconds  x (number  of  maps) 


C+'J. 
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SECTION  VI 
PROGRAM  LIMITATIONS 


Program  PHASE4  requires  that  the  number  of  vehicle-identification  cards 
be  from  one  to  10  because  of  the  array  dimensions.  The  number  of  map-bounds 
cards  may  range  from  0 to  100,  but  the  number  of  sections  without  bounds  cards 
plus  the  number  of  bounds  cards  must  not  exceed  100.  Array  dimensions  impose 
this  limitation,  also. 

The  plotter  may  be  either  a flatbed  or  a drum  plotter  with  a height  of  at 
least  30  inches.  If  a smaller  drum  plotter  is  used,  the  current  assignment 
YHCUT=30. , which  appears  shortly  after  statement  330  in  the  main  program,  should 
be  replaced  by  a statement  assigning  to  YHCUT  the  value  of  the  drum  height,  in 
inches. 

The  use  of  a colon,  a nonstandard  CDC  character,  in  the  printing  of  time 
imposes  a restriction  on  the  storage  and  transmittal  of  the  FORTRAN  listing 
of  the  program  and  the  route  schedule.  The  colon  will  survive  as  a colon  only 
on  disk  or  cards.  If  a route  schedule  or  the  card  images  of  the  program  are 
transferred  to  magnetic  tape,  the  colon  will  be  converted  to  a zero.  On  some 
devices,  such  as  a microfiche  unit,  the  colon  may  be  completely  deleted.  One 
colon  is  used  on  each  of  the  cards  with  serial  numbers  PHS42700,  HM001100,  and 
PRSC0220.  It  will  be  necessary  to  reinsert  the  colon  on  these  cards  if  the 
cards  are  transferred  to  tape. 


I 
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SECTION  VII 

ERROR  MESSAGES  AND  CORRECTIVE  ACTION 


1.  NO  MAP  BOUNDS  WERE  GIVEN  FOR  ANY  TRIP. 

EACH  MAP  WILL  SHOW  TRAVEL  IN  THE  COLLECTION  REGION  BUT  NOT 
NECESSARILY  THE  PATH  TO  OR  FROM  THE  GARAGE  OR  LANDFILL. 

Type:  Warning. 

Source:  Program  PHASE4. 

Location:  Where  map-bounds  listing  would  normally  appear. 

Meaning:  Self-explanatory. 

Action:  If  maps  including  travel  to  or  from  the  garage  or  landfill  are 

desired,  use  map-bounds  cards. 


2.  WHERE  NO  BOUNDS  WERE  SPECIFIED,  THE  MAP  WILL  SHOW  TRAVEL  IN  THE 
COLLECTION  REGION  BUT  NOT  NECESSARILY  THE  PATH  TO  OR  FROM  THE 
GARAGE  OR  LANDFILL. 


Type: 

Source: 

Location: 

Meaning: 

Action: 


Informative. 

Program  PHASE4. 

At  the  end  of  the  map-bounds  listing. 

Self-explanatory. 

If  travel  to  or  from  the  garage  or  landfill  is  desired  on  all  maps, 
use  map-bounds  cards. 


3.  NO  PATH  DATA  FOUND  ON  UNIT  9 OR  ON  CARDS. 

JOB  TERMINATED. 

Type:  Fatal. 

Source:  Program  PHASE4. 

Location:  After  map-bounds  listing. 

Meaning:  No  path  data  were  found  on  cards  or  on  file  TAPE9. 

Action:  If  the  path  data  should  be  on  TAPE9,  see  whether  an  ATTACH  or 

REQUEST  card  for  that  file  is  present  and  correct.  If  path  data 
should  be  on  cards,  look  for  an  extra  end-of-record  card  or  an 
end-of-file  card  before  the  data.  The  deck  may  be  missing. 
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4.  THE  MAXIMUM  TRIP  TIME  WILL  BE  EXTENDED  TO  hh.hh  HOURS  FOR 
VEHICLES  OF  CAPACITY  vvvvv.v 

IF  THIS  IS  UNSATISFACTORY,  RERUN  PROGRAM  PHASE2  and  PHASE3 
WITH  A SMALLER  TIME  LIMIT  IN  PHASE2. 

Note:  hh.hh  represents  a time  in  decimal  hours, 

vvvvv.v  represents  a vehicle  capacity. 

Type:  Warning. 

Source:  Program  PHASE4. 

Location:  Following  the  printing  of  the  trip  descriptions. 

Meaning:  A vehicle  with  the  indicated  capacity  requires  more  than  the 

user-specified  maximum  trip  time.  The  program  will  automatically 
extend  the  maximum  trip  time  to  allow  longer  trips. 

Action:  If  the  time  extension  is  not  acceptable,  rerun  programs  PHASE2 

and  PHASE3,  using  a smaller  time  limit  in  PHASE2. 


Meaning: 


Action: 


5.  PIECE  p,  nnnTH  NODE,  NUMBER  mmmmm,  IS  INCORRECT. 

Note:  p is  either  1,  2,  or  3. 

nnn  is  a sequence  number  from  1 through  101, 
mmmmm  is  a node  number  from  1 through  99999. 

Type:  Warning,  but  the  schedule  will  be  incorrect. 

Source:  Subroutine  PATHPLT. 

Location:  Following  the  printing  of  the  schedule. 

Meaning:  Node  number  mmmmm,  used  in  the  path  description  on  cards  or  on 

file  TAPE9,  is  not  present  in  the  node  data  from  file  TAPE2. 

Action:  Verify  that  file  TAPE2  and  the  path  data  on  TAPE9  or  cards  corre- 

spond to  the  same  problem.  If  the  path  is  on  cards,  check  for  a 
card-punch  error.  If  the  total  number  of  steps  for  the  three 
pieces  exceeds  200,  the  problem  originates  in  program  PHASE3.  The 
last  two  paragraphs  in  Volume  III,  Section  VII,  of  this  report  de- 
scribe the  cause  of  the  error  and  possible  corrective  actions. 

6.  PIECE  p,  nnnTH  SEGMENT,  NUMBER  mmmm,  DOES  NOT  CONNECT  TO  A 
BOUNDING  NODE. 

Note:  p is  either  1,  2,  or  3. 

nnn  is  a sequence  number  from  1 through  100. 
mmmm  is  a segment  number  from  1 through  1023. 
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Type:  Warning,  but  the  schedule  will  be  incorrect. 

Source:  Subroutine  PATHPLT. 

Location:  Following  the  printing  of  the  schedule. 

Meaning:  One  or  both  nodes  bounding  segment  imrnn  in  a path  description 

do  not  bound  the  segment  in  the  data  from  file  TAPE1. 

Action:  Verify  that  file  TAPE1  and  the  path  data  on  TAPE9  or  cards  corre- 

spond to  the  same  problem.  If  the  path  is  on  cards,  check  for  a 


card-punch  error.  If  the  total  number  of  steps  for  the  three 
pieces  exceeds  200,  the  problem  originates  in  program  PHASE3.  The 
last  two  paragraphs  in  Volume  III,  Section  VII,  of  this  report  de- 
scribe the  cause  of  the  error  and  possible  corrective  actions. 

7.  NO  VEHICLE  SPECIFIED  FOR  THIS  CAPACITY. 

Type:  Warning,  but  the  schedule  may  be  incorrect. 

Source:  Main  program  PHASE4. 

Location:  In  the  route  summary  and  in  the  heading  for  the  schedule. 

Meaning:  A vehicle  capacity  from  the  route  data  on  either  TAPE9  or  cards 

is  not  included  in  a vehicle-identification  card. 

Action:  Check  the  vehicle-identification  cards.  A vehicle  capacity  may 

be  missing  or  mispunched. 
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SECTION  VIII 

RECOMMENDED  PROGRAM  CHANGES 


Two  changes  in  program  PHASE4  are  recommended.  The  trip-pairing  algorithm 
in  the  loop  through  statement  620  in  the  main  program  could  be  modified  to  gen- 
erate single-trip  routes  when  a two-trip  route  would  exceed  the  maximum  working 
time  per  day  (TMXDAY)  specified  by  the  user.  The  change  would  require  moderate 
modifications  of  the  coding  within  the  loop,  and  may  require  minor  modifications 
in  the  path-,  map-,  and  schedule-generation  subroutines. 

At  present,  the  sides  of  the  streets  are  plotted  shorter  than  the  segment 
length  so  that  they  meet  only  where  streets  come  together  at  right  angles.  The 
fourth  and  fifth  arguments  to  subroutine  SHAPCOM  allow  length  modifications  at 
the  initial  and  final  ends  of  the  segment.  Additional  coding  could  be  inserted 
before  the  call  to  SHAPCOM  in  subroutine  MAPPLT  to  adjust  street  lengths  to  take 
into  account  the  angle  between  the  segments.  The  change  would  require  a know- 
ledge of  geometry,  trigonometry,  and,  possibly,  vectors. 
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Subroutine  STRIN 
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Convert  whole  and  fractional 
hours  to  hours  and  minutes, 
rounding  to  the  nearest  minute. 
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Set  initial  and  final  pointers. 


Save  values  at  inspection  point. 
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Subroutine  SHLSRT 
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Subroutine  P0SIT9 
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Set  straight  segment  parameters. 


Compute  slope  and  distance 
between  segment  endpoints. 


Compute  parameters. 
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Print  travel  or 
collection  instructions. 


Print  break  instructions. 


Print  lunch  instructions. 
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Call  SHLSRT  to  sort  the  segments 
by  length,  the  longest  first. 


Subroutine  STNAME 
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END 


Call  COORD  and  PLOT  to 
plot  path  (solid  line 
indicates  travel,  only; 
dashed  line  indicates 
collection). 


Program  PHASE4 
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Call  CUMTD  to 
accumulate 
time  and  distance 
of  trip  from 
landfill  to  garage. 


Read  trip-path  description  data  from  TAPE9  or  cards 
(TAPE5),  as  determined  previously,  for  trips  to 
region,  within  region,  and  from  region. 


If  trip  description  is 
read  from  cards,  write 
data  to  TAPE9. 


Call  CUMTD  to 
accumulate  times 
and  distances  of 
trio  pieces. 
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Call  SHLSRT  to  sort  trips 
by  vehicle  capacity. 


Call  SHLSRT  to  sort  trips 
having  the  same  vehicle 
capacity  by  time. 
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Pair  an  afternoon  and  a 
morning  trip  to  make  a 
route. 


Generate  totals  to  be  used 
in  printing  a route  summary. 


Program  PHASE4 
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Function  I FIND 
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Subroutine  STRIN 
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Function  HM 
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Subroutine  NUMBER 
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Subroutine  CUMTD 
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Subroutine  SHLSRT 

131 

Subroutine  P0SIT9 

131 

Subroutine  SHAPCOM 

131 

Subroutine  COORD 
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Subroutine  PRSCHED 

132 

Subroutine  STNAME 

133 

Subroutine  MAPPLT 
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Subroutine  PATHPLT 

135 

Program  PHASE4 

136 

Note:  A single  variable  symbol  may  have  different  meanings  in  relation  to  the 
various  subroutines.  For  this  reason,  variables  are  defined  below  for 
each  subprogram  and  for  the  main  program. 


I ARRAY 


FUNCTION  IFIND 


Array  being  searched. 
Length  of  IARRAY. 
Number  being  sought. 


NAMSTR 

NSTRS 

NUMSTR 


SUBROUTINE  STRIN 

Array  of  street  names. 

Count  of  street-name  cards. 

Array  of  street  numbers. 

FUNCTION  HM 

Integer  number  of  hours  in  TIME. 
Number  of  minutes. 

Time,  in  decimal  hours. 


SUBROUTINE  NUMBER 

Output  format  for  number. 

Number  to  be  plotted. 

Character  representation  of  number. 


SUBROUTINE  CUMTD 

Array  of  collection-or-travel  indicators. 

Length  of  path,  in  miles. 

Array  of  segment  lengths,  in  miles. 

Array  of  segment  speed  limits. 

Array  of  numbers  of  segments  in  path. 

Array  of  segment  house  counts. 

Total  house  count  for  path. 

Count  of  segments  in  path. 

Total  refuse  quantity  for  path. 

Array  of  segment  refuse-quantity  adjustment  factors. 


130 


(mm 


SCOLL  Vehicle  speed  during  collection. 

TIM  Time  required  to  traverse  path,  in  minutes. 

TSTOPH  Stop  time  per  collection  location,  in  minutes. 

TSTOPR  Stop  time  per  unit  refuse,  in  minutes. 


A 

NW 

SGN 

X 


SUBROUTINE  SHLSRT 

Array  reordered  as  X is  sorted. 

Count  of  words  to  be  sorted. 

If  SGN  = 1.0,  X is  sorted  into  increasing  order. 
If  SGN  = -1.0,  X is  sorted  into  decreasing  order. 

Array  to  be  sorted. 


KLIM 

KS 

KSKP 

LSEQN 

NSEQN 


SUBROUTINE  P0SIT9 

Count  of  pieces  in  trip. 

Count  of  segments  in  piece. 

Count  of  card  images  to  be  skipped. 

Sequence  number  of  last  trip  read  from  TAPE9. 
Sequence  number  of  trip  being  sought  on  TAPE9. 

SUBROUTINE  SHAPCOM 


- 


A 

AVMD 

AW 

B 

BR1 

BR2 

DF 

DI 

DIR 

DX 

DY 

ISF 

RATIO 


. _ . _ 


Array  giving  slopes  of  linear  mappings  of  segment  pieces. 

Map  distance  conversion,  in  miles  per  MCU. 

Absolute  value  of  street  half-width,  in  miles. 

Array  of  constants  in  the  linear  mappings  of  segment  pieces. 

Distance  to  first  break  in  segment  shape,  in  miles. 

Distance  to  second  break  in  segment  shape,  in  miles. 

Displacement  of  end  point  of  segment  mapping,  in  miles. 

Displacement  of  starting  point  of  segment  mapping,  in  miles. 

If  DIR  =1.0,  coordinates  are  generated  from  the  starting  to  the 
ending  node.  If  DIR  = -1.0,  coordinates  are  generated  from  the 
ending  to  the  starting  node. 

X-component  of  segment  span,  in  MCU. 

Y-component  of  segment  span,  in  MCU. 

Shape  code  when  in  character  form. 

Array  of  ratios  of  side  length  to  segment  length  for  segment 
pieces. 
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a .. 


) 


| 

«■ 


RPR 

SF 

THETA 

TOTLEN 

W 

XNF.YNF 

XNI.YNI 


Reciprocal  of  radius  of  curvature  for  circular  segments  or 
circular  portions  of  S-curves. 

Shape  code  when  in  binary  form. 

Slope  of  line  from  starting  to  ending  node,  in  radians. 
Segment  length,  in  miles. 

Half-width  of  segment,  in  miles. 

Coordinates  of  ending  node. 

Coordinates  of  starting  node. 


A 

B 

BR1 

BR2 

CUMLEN 

NPC 

RATIO 

RPR 


S 

SF 

XINT.YINT 


SUBROUTINE  COORD 

Array  of  slopes  of  linear  mappings  of  segment  pieces. 

Array  of  constants  in  the  linear  mappings  of  segment  pieces. 
Distance  to  first  break  in  segment  shape,  in  miles. 

Distance  to  second  break  in  segment  shape,  in  miles. 

Cumulative  length  along  segment,  in  miles. 

Number  of  piece  of  segment. 

Array  of  ratios  of  side  length  to  segment  length  for  segment 
pieces. 

Reciprocal  of  radius  of  curvature  for  circular  segments  or 
circular  portions  of  S-curves. 

Linearly  mapped  distance  along  segment  since  previous  break. 
Shape  code. 

Arrays  of  coordinates  used  in  linear  interpolations  for  segment 
coordinates. 


XNF.YNF 

XNI.YNI 


Coordinates  of  ending  node. 
Coordinates  of  starting  node. 


ACT 

CORT 

DBRK 

DIS 

DLUNCH 

FLEN 

FMPH 

KNODES 


SUBROUTINE  PRSCHED 

Collection-or- travel  indicator  for  current  segment. 

Array  of  collection-or- travel  indicators  for  the  segments  in 
the  path. 

Array  of  durations  of  break  times,  in  minutes. 

Distance  traveled  since  last  line  of  printed  output,  in  miles. 
Duration  of  lunch  break,  in  minutes. 

Array  of  segment  lengths,  in  miles. 

Array  of  segment  speed  limits,  in  mph. 

Count  of  nodes  in  map  description. 
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MXLINE 

NAMSTR 

NFS 

NH 

NMSG1 

NMSG2 

NNP 

NSG 

NSP 

NSTR 

NSTRS 

NTPS 

NUMSTR 

OLDTIM 

RQF 

TBRK 

TIM 

TLUNCH 

TOTLD 

TSTART 

TSTOPH 

TSTOPR 

TUNLD 


AHGT 

AVMD 

AWDTH 

CORT 

FLTEM 

ICHAR 

INS 

KNODES 

NAMSTR 

NCH 


Count  of  lines  per  page  before  page  eject. 

Array  of  street  names. 

Count  of  stops  since  previous  line  of  printed  output. 

Array  of  segment  house  counts. 

Number  of  current  segment. 

Number  of  next  segment  in  path. 

Array  of  numbers  of  nodes  in  path. 

Array  of  counts  of  segments  in  path  pieces. 

Array  of  numbers  of  segments  in  path. 

Array  of  street  numbers  of  segments. 

Count  of  street  names. 

Count  of  trip  choices  available  per  section. 

Array  of  street  numbers. 

Time  at  completion  of  previous  segment. 

Array  of  segment  refuse-quantity  adjustment  factors. 

Array  of  approximate  break  starting  times. 

Time  at  end  of  current  segment. 

Approximate  starting  time  of  lunch  break. 

Cumulative  vehicle  load. 

Route  starting  time,  in  hours. 

Stop  time  per  collection  location,  in  minutes. 

Stop  time  per  unit  of  refuse,  in  minutes. 

Unloading  time  at  the  landfill,  in  minutes. 

SUBROUTINE  STNAME 

Height  of  the  lettering,  in  inches. 

Map  distance  conversion,  in  miles  per  MCU. 

Width  of  the  street  name,  in  inches. 

Array  of  collection-or-travel  indicators  for  segments  in  path. 
Array  of  segment  lengths,  in  miles. 

Array  of  characters  in  street  name. 

If  INS  = 0,  street  name  is  plotted  along  outside  of  street; 
if  INS  = 1,  street  name  is  plotted  within  street  bounds. 

Count  of  nodes  in  map  description. 

Array  of  street  names. 

Count  of  characters  in  street  name. 
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NMAP 

NMAPO 

NN1 

NN2 

NODNUM 

NSEG 

NSG 

NSP 

WIDTH 

XL  ,XR 

XNF.YNF 

XNI.YNI 

XNOD.YNOD 

XPF.YPF 

XPI.YPI 

YB.YT 

YHCUT 


AVMD 

CUMLEN 

FLEN 

INB 

ISF 

KNODES 

NMAP 

NMAPO 

NN1 

NN2 

NODNUM 

NPPSEG 

PHGT 

PLEN 

SF 

TOTLEN 

W 


T 


Map-strip  number  of  the  current  point. 

Map-strip  number  of  the  previous  point. 

Array  of  segment  starting  node  numbers. 

Array  of  segment  ending  node  numbers. 

Array  of  node  numbers. 

Count  of  segments  in  map  description. 

Array  of  counts  of  segments  in  each  piece  of  the  trip. 

Array  of  numbers  of  segments  in  trip. 

Half-width  of  the  street,  in  inches. 

X-coordinates  of  left  and  right  map  boundaries,  in  MCU. 
Coordinates  of  ending  node,  in  MCU. 

Coordinates  of  starting  node,  in  MCU. 

Arrays  of  node  coordinates. 

Plotter  coordinates  of  start  of  next  character  in  street  name. 
Plotter  coordinates  of  start  of  present  character  in  street  name. 
Y-coordinates  of  bottom  and  top  map  boundaries,  in  MCU. 

Height  of  map  strips,  in  inches. 


SUBROUTINE  MAPPLT 

Map  distance  conversion,  in  miles  per  MCU. 
Cumulative  length  along  segment,  in  miles. 
Array  of  segment  lengths,  in  miles. 

Point  within  map-bounds  indicator. 

Shape  code  when  in  character  form. 

Count  of  nodes  in  map  description. 
Map-strip  number  of  the  current  point. 
Map-strip  number  of  the  previous  point. 
Array  of  starting  node  numbers. 

Array  of  ending  node  numbers. 

Array  of  node  numbers. 

Count  of  points  plotted  per  segment. 

Height  of  map  strip,  in  inches. 

Total  length  of  plot,  in  inches. 

Shape  code  when  in  binary  form. 

Segment  length.  In  miles. 

Half-width  of  street,  in  miles. 
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T 


WIDTH 
XL  ,XR 
XLEN 

XNOD.YNOD 

YB.YT 

YCUT 

YLEN 


ACT 

AVMD 

CUMLEN 

FLEN 

INB 

ISEG 

ITRV 

KARO 

KK 

NMAP 

NMAPO 

NNP 

NPPSEG 

NSG 

NSP 

PHGT 

RSO 

SF 

TOTLEN 

W 

WIDTH 

XL,XR 

XLEN, YLEN 

XNOD.YNOD 

YB.YT 

YCUT 


Half-width  of  street,  in  plotter  inches. 

X-coordinates  of  left  and  right  map  boundaries,  in  MCU. 
X-direction  map  length,  in  inches. 

Arrays  of  node  coordinates,  in  MCU. 

Y-coordinates  of  bottom  and  top  map  boundaries,  in  MCU. 
Height  of  map  strips,  in  MCU. 

Y-direction  map  length,  in  inches. 

SUBROUTINE  PATHPLT 

Col lection-or-travel  indicator  for  current  segment. 

Map  distance  conversion,  in  miles  per  MCU. 

Cumulative  length  along  segment,  in  miles. 

Array  of  segment  lengths,  in  miles. 

Segment-in-bounds  indicator. 

Array  of  numbers  of  segments  in  path. 

Array  of  counts  of  times  segments  are  traversed. 

Count  of  steps  before  an  arrowhead  is  plotted. 

Number  of  current  segment. 

Map-strip  number  of  current  point. 

Map-strip  number  of  previous  point. 

Array  of  numbers  of  nodes  in  path. 

Count  of  points  plotted  per  segment. 

Array  of  counts  of  segments  in  pieces  of  trip. 

Array  of  numbers  of  segments  in  path. 

Height  of  map  strip,  in  inches. 

Indicator  for  collection  from  only  right  side. 

Shape  code. 

Segment  length,  in  miles. 

Half-width  of  street,  in  miles. 

Half-width  of  street,  in  plotter  inches. 

X-coordinates  of  left  and  right  map  boundaries,  in  MCU. 
X-  and  y-direction  map  lengths,  in  inches. 

Arrays  of  node  coordinates,  in  MCU. 

Y-coordinates  of  bottom  and  top  map  boundaries,  in  MCU. 
Height  of  map  strips,  in  MCU. 
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PROGRAM  PHASE4 


AVMD 

CORT 

DBRK 

DIST 

DLUNCH 

DTIF 

INPTU 

IORD 

IPAIR 

IRS 

KNODES 

LV 

NBC 

NHS 

NNP 

NRT 

NSC 

NTC 

NTPS 

NRTP 

SC 

SCOLL 

TBRK 

TC 

TCAP 

TLOAO 

TLUNCH 

TMXDAY 

TMXTR 

TSTART 

TSTOPH 

TSTOPR 

TTIF 


Map  distance  conversion,  in  miles  per  MCU.  ' 

Array  of  col lection-or-travel  indicators.  \ 

Array  of  break  durations,  in  minutes. 

Array  of  distances  to  section,  within  section,  to  landfill  from 
section,  and  from  landfill  to  garage. 

Duration  of  lunch  break,  in  minutes. 

Array  of  travel  distances  to  section,  within  section,  and  from 
section  to  landfill,  for  all  sections. 

Number  of  input  unit  containing  route  data. 

Array  of  pointers  to  trips. 

Consecutive  trip-pairing  indicator. 

Array  of  section  numbers  in  order  of  output. 

Count  of  nodes  in  map  description. 

Array  of  vehicle-identification  line  numbers. 

Count  of  map-bounds  cards. 

Array  of  house  counts  per  section. 

Array  of  numbers  of  nodes  in  path. 

Route  counter. 

Section  number. 

Count  of  vehicle  identifications. 

Count  of  trip  choices  per  section. 

Array  of  trip  indicators  from  map-bounds  cards. 

Map  scale,  in  inches  per  MCU. 

Vehicle  speed  during  collection. 

Array  of  approximate  break  starting  times,  in  hours. 

Array  of  vehicle  capacities  from  vehicle-identification  cards. 

Array  of  capacities  of  vehicles  servicing  sections. 

Array  of  vehicle  loads. 

Approximate  starting  time  of  lunch  break,  in  hours. 

User-specified  maximum  route  time  per  day. 

User-specified  maximum  trip  time,  in  hours. 

Route  starting  time,  in  hours. 

Stop  time  per  collection  location,  in  minutes. 

Stop  time  per  unit  of  refuse,  in  minutes. 

Array  of  travel  times  to  section,  within  section,  and  from 
section  to  landfill,  for  all  sections. 

|Lj*^.3k 
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TUNLD  Unloading  time  at  the  landfill,  in  minutes, 

UNITS  Refuse  unit. 

VID  Array  of  vehicle  identifications. 

WIDTH  Half-width  of  street,  in  plotter  inches. 

XNOD.YNOD  Arrays  of  node  coordinates,  in  MCU. 
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APPENDIX  E 

SAMPLE  ROUTE  MAPS 
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APPENDIX  F 

SAMPLE  PRINTED  OUTPUT 
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GLOSSARY 

Air  Force  Refuse-Collection  Scheduling  Program:  a set  of  four  computer  pro- 
grams that  schedule  residential  refuse-collection  and  produce  printed 
schedules  and  maps  of  the  routes. 

map  coordinate  unit  (MCU):  the  length,  in  inches,  of  a unit  interval  from  the 
coordinate  system  appended  by  the  user  to  the  first  map  input  to  pro- 
gram RCINPT. 

node:  a numbered  point  on  a street  at  which  some  characteristic  of  the  street 
changes. 

route:  the  travel  and  collection  performed  by  one  collection  vehicle  during  a 
day,  starting  and  ending  at  the  garage. 

section:  a group  of  segments  serviced  by  the  same  collection-vehicle  trip. 

segment:  a portion  of  a street  between  two  nodes. 
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